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A number of microcrystal layers of evaporated lead sulphide in vacuum envelopes 
were investigated. The previously developed method (H. Rzewuski, L. Sosnowski 1955) 
of investigating semi-conductor layers by applying, in a condenser-like arrangement, 
a pulsed electric field to the layer is discussed extensively. On the basis of the above method, 


the following were determined: 


a. 


5. 


ово г 


the 


. the 


the 
the 
the 
the 
the 


type of layer investigated 

type of barrier layer 

concentration of the majority carriers 

mobility of the majority carriers 

minimum density of the surface states 

relaxation time of the barrier 

kinetic changes of the barrier layer conductivity are discussed. 


Finally the accuracy of the pulse method used to investigate the type of semi-con- 
ductor layer, as well as possibilities of its application, is discussed. 


1. Surface States 


Part I 


The intensive development of the field of semi-conductors which has taken place 
particularly in the last few years is connected with the discovery of many interesting 
properties of semi-conductors having a basic importance to electronics. This had 
lead, in turn, to the extension of the present theoretical picture of semi-conductors 
and has carried the main weight of the investigations to non-homogeneous arrange- 
ments with strongly variable concentrations of current carriers. A role of prime impor- 
tance is played here by {һе contact problem, the related question of the n-p barrier, 


and surface problems. 


From the great deal of work done in this field, it has been established that the 
contact phenomena basically depend on the character of the surfaces of the contacting 
_ bodies and on the density of the surface states and their energy distribution. In reality, 
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; 


958 H. Rzewuski 


the picture.of the phenomena taking place on the surface of a semi-conductor is very 
complex. Generally speaking, wherever the phenomena have a surface character, the 
surface states play the decisive role. 

In the band model of the solid body we have to do with a simplifying assumption, 
in which we consider a crystal to be boundlessly extended along the three axes x, y, z 
In this concept, we ‘lose‘ the entire class of phenomena connected with the crystal 
surface, which constitutes a disturbance of the periodicity of the lattice by the very 
fact of its being bounded. Proper conclusions from this fact were drawn by Tamm 
(1932, 1955) who obtained, with a simplified crystal model (the Kronig-Penney model), 
a one-dimensional — and subsequently three-dimensional — group of electron states 
localized on the surface and lying in the energy interval between the basic band and 
the conductivity band. After some work in this field (S. Ryzhanow 1934, A. Sokolov 
1934, 1936, A. Maue 1935, E. Goodwin 1939, R. Fowler 1953) a major step forward 
in the theory of surface states was made by Shockley (1939). He examined the change 
in the wave functions of the electron and the electron states with the change in intera- 
tomic distances of a one-dimensional crystal lattice composed of a small number of 
atoms. For large lattice constants the energy levels are the same as for the individual 
atoms, but for smaller distances we obtain a distinct degeneration of the energy levels 
found in the allowed energy band. With a further decrease of the lattice constant 
the allowed bands begin to overlap and there appear in the forbidden energy regions 
two levels corresponding to the local surface states. Extending these considerations 
to three dimensions, we obtain in the energy gap a band of surface states. 

At present the notion of surface states is understood in a sense broader than the 
normal Tamm states. The expression ‘disturbance of the periodicity of the lattice‘ 
at the surface is extended also to all disturbances, both mechanical (dislocations, detects) 
and chemical (acceptor and donor centres) occurring in the surface layer. The thickness 
of this layer is taken to be of the order of 1077 cm, i.e. the free path of an electron 
(because of the treatment received by the surface and the action of the surrounding 
medium, the surface layer may have properties other than those inside the crystal). 
The type of potential barrier formed just under the semi-conductor surface will depend 
on the character of the surface states and on their density. For a large enough density 
of surface states, there will arise on the free surface of the semi-conductor a double 
layer of charge (1947) consisting of: 5 


a. a surface layer (of the order of 10-7 ст) ете from electrons (or holes) . 


localized in surface states; 

b. a layer of space charge of opposite sign extending from the surface of the 

crystal to a depth of 10-9 to 10-* cm. 

А non-homogeneous charge distribution in the vicinity of the surface and connec- 
ted with the formation of a double layer, causes the electron levels to overlap in the 
space charge region and creates the so-called potential barrier. 

This paper is based on the phenomenon of the injection of barriers by an electric 
field applied perpendicularly to the surface of the investigated layer in a condenser 
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arrangement. The effect of the electric field on the conductivity of the semi-conductor 
was observed for the first time by Shockley and Pearson in 1948 in connexion with 
the investigation of germanium. It turned out, however, that the modulation of the 
conductivity of the semi-conductor is much smaller than expected. Bardeen (1947) 
showed that a considerable part of the induced carriers are localized on the surface 
and pass over to surface states. In this work he showed the reason for the lack of 
a correlation between the contact potential and the difference of work functions of the 
metal and the semi-conductor as observed by Meyerhoff (1947). The presence of surface 
states influences the overlapping of the energy bands at the surface; in the case of a large 
surface state density a double layer forms on the surface of the semi-conductor from 
the charge in the surface states and the charge of opposite sign in the barrier layer 
just under the surface. The potential difference where the metal is in contact with the 
semi-conductor is compensated mainly by the surface charge, and аз а result, the barrier 
potential does not undergo change. The existence of a surface conductivity resulting 
from the charges in the barrier layer, this conductivity being of a sign opposite to the 
volume conductivity, was confirmed in the work of Brown (1953) in which.he showed 
that an n-type surface layer forms on a sample of p-type germanium. 


2. Method of Measurement 


А diagram of the arrangement used by Shockley and Pearson to modulate the 
conductivity of the germanium layer by means of surface charges given below. 

Direct current flows in the circuit W К U which includes the layer under investi- 
gation (Fig, 1). The value of this current depends, in part, on the resistance of the 
layer W. Application of an electric field perpendicularly to the surface of the layer 
in a condenser arrangement, where the metal electrode forms one plate and the layer 
itself, the other, results in an induced charge q, on the surface of the layer; this positive 
charge changes the conductivity of the entire layer. The change of conductivity may 
be recorded by observing either the change in current in the circuit W R U, or the 
change of voltage across resistance .R. Because of the existence of electronic surface 
_ states on the surface of the layer, part of the charge g, is immobilized іп the surface 
—states, and only the difference q, = qo — q, enters into the barrier layer lying just 
below the surface layer. The electrons localized on the surface do not take part in the 
conductivity, while the electrons introduced into the barrier layer pass over the con- 
ductivity band and change the local concentration of the carriers, the conductivity 
of the entire layer now being changed to a lesser degree. It is precisely this change 
of conductivity which is recorded and which forms the basis of the determination of 
the density of the surface states. It should be stressed that the above method is based 
on the observation of stationary- steady-state effects. 


In the method worked out, it was attempted to obtain the largest possible amount 
of information concerning the semi-conductor under investigation without resorting 


i 
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to supplementary measurements.: The characteristic features of the method, drawn 
up together with Professor L. Sosnowski (1955), are: 

a. the applied electric field had the form of rectangular pulses 

b. a bridge measuring circuit allowed compensation for spurious signals 

c. in connexion with the above, it was possible to investigate the effect for the 

non-steady state case. 

The semi-conductor elements investigated in this work consisted of a layer of 

microcrystalline PbS evaporated in, vacuo on glass between two very thin graphite 


Tungsten contacts 


03-08 


Fig. 1 | Fig. 2 


electrodes. The method of producing the layers was based on the work of Professor 
L. Sosnowski (1949). The approximate dimensions of the layers and the glass vacuum 
envelope in which they were evaporated are given below (Fig. 2). 

The active element here is the layer contained between the graphite electrodes. 
The electric field is applied in a condenser arrangement through the walls of the envelo- 
pe. One plate is formed by the metal electrode pressed against the bottom of the 
envelope, while the second is the layer itself. 


Investigation of the above layers by means of an alternating current field in the 
Shockely-Pearson arrangement was not possible because of the unavoidable capacitive 
coupling between the metal electrode and the graphite electrodes. Approximate cal- 
culations already show that the spurious signal of capacitive origin not only makes 
the measurement of the effect itself difficult but renders it completely impossible. 
In the case of the investigated layers, the stray capacitance had the value of several pF, | 
which, in conjunction with the high resistance of the layer — hundreds of kQ (and 
therefore also the matching resistance А) — gave a spurious signal for the lowest 
possible frequencies of the applied voltages, this signal exceeding the expected effect 
by about three orders of magnitude. The situation was made difficult by the fact that 
the applied field had a rectangular wave shape, and therefore contained a relatively 
broad frequency spectrum of various amplitudes (Fourier distribution). Thus the com- 
pensation used had to contain amplitude and phase corrections over the entire trans- 
mitted frequency band without distorting the measured effect. 


A schematic diagram of the measuring arrangement is shown below іп Fig. 3. 
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1. Rectangular wave generator with negative polarity gives a voltage of U п. = 
= 0—100 V. The pulse repetition rate can be changed from 5 to 50000 c/s and the 
pulse width from 19/20 to 1/20 of a cycle. The rise and fall time of the pulse is not 
worse than 0.5 usec. 

II. Cell with the investigated layer works in a bridge circuit with adjustable pulse 
amplitude and phase. 

Two circuits should be distinguished here (Fig. 3): one circuit, working with 
the cell includes the layer resistance R,, the adjustable resistances R, and R,, and 
the battery U_; in the second circuit the generator voltage is fed to the stray capacities 
Cp, and С, (shown in Fig. 4), variable capacities C, and Cj, resistors R, and R,, and 
battery U_ to earth. The procedure is as follows. The voltage (п. from the generator 
is continuously connected to the metal electrode. At the beginning U_ is disconnected 
from the circuit (the switch is not shown in the schematic diagram) and, as a result, 


Fig. 3 Fig. 4 


the effect cannot occur. At outputs A, and A, we have only the spurios signal. By 
regulating the resistors R, and R, and the capacities C, and C, we adjust the phase 
and amplitude of the bridge so as to obtain a minimum signal at the output (i.e. on 
the oscillograph screen). We now connect the battery U_ to the cell circuit and, because 
of the pulsed field, we also obtain a voltage derived from the effect under investigation. ' 
At the outputs A, and A,, symmetric with respect to earth, we have the signals resulting 
from the change of conductivity — of the order of mV — and spurious signals — of ` 
the order of volts. The voltage is fed to a differential circuit with a symmetric cathode 
follower input and then to a symmetric oscillograph of the *Furzehill" type. The last 
two stages, i.e. III and IV (differential amplifier and oscillograph) react only to the 
difference of the signals arriving at their inputs. Signals of the same amplitude and 
phase with respect to earth are not recorded by the final circuit. In the case discussed, 
the spurious signal did not cause any deflexion of the oscillograph when the bridge 
was carefully balanced, since at both inputs 4, and A, we have voltages of precisely 
the same amplitude and phase. On the other hand, the changes in conductivity of the 
layer produce at inputs A, and A, voltages of the same amplitude, but of opposite 
phases with respect to earth. As a result, the voltages on the oscillograph plates add 
to each other and we observe a full deflexion on the screen. Resistances R, and R, 
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can be adjusted so that their values are equal while their sum constitutes a matching 
resistance to the layer,.i.e. R, + R, = Rẹ The bridge elements are so chosen that 
the bridge may be rapidly adjusted. Any eventual inequality which may arise in 
adjusting resistances R, and R, simultaneously does not have any essential influence 
on the resistance matching of the cell resistance since the curve is sufficiently flat. 

Ш. This part of the arrangement constitutes a two stage symmetric circuit with 
an input cathode follower using 6AK5 pentodes operating as triodes. А high resistance 
input is required because of the large resistance 
of the layers. As is known, the dynamic resis- 
tance of such a circuit is approximately 


Rsk 


Rayn = ЖЕ Ой” (1) 


where R,, is the grid-to-cathode resistance, and 
k„ is the amplification of the valve in the ca- 
thode follower circuit. The schematic diagram of 
this stage is shown in Fig. 5. 

The next stage is a discriminator circuit with 
a small amplification. The coefficient of discrimi- 
nation, taken here as the ratio of the output to 
the in-phase input signal, is high. 

IV. The voltage taken symmetrically with 
the output of this circuit, is next fed 10 
the oscillograph, which also has a high coefficient of discrimination and a large 
amplification of signals of opposite phase. The circuit as a whole permits the measure- 
ment of effects of the order of millivolts. As a result of the compensation the spurious 
signals are reduced to a level several times lower than the measured effect. 


4. Oscillograph Wave Forms 


We shall now examine the oscillograph wave forms of the conductivity changes 
in the layer as a result of the application of the pulse field. (By the term ‘pulse field‘ 
we will henceforth understand the field consisting of rectangular pulses.) Let the pulse 
under consideration be sufficiently far in time from the instant the circuit is switched in. 

At 1 (Fig. 6) the conductivity of the layer is smaller than for the case of no field. 
The conductivity could not attain the value for the case of a layer without a field, since 
the following impulse arrived. Application of the pulse field at a time corresponding 
to 1 causes an increase in charge in the layer volume and thus an increase in the volume 
conductivity (which changes with changes in the field). The instantaneous conductivity 
increases to the value corresponding to 2 (it should be noted that the wave forms obtain- 
ed correspond to changes in the conductivity and not to its total value). The increase 
in conductivity of the layer in relation to the value without the field — dotted line — 
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corresponds to the segment 1’—2, since part of the change in conductivity went to 
compensate the value 1—1’. The pulse field now maintains a constant amplitude, but 
the increase in conductivity slowly reduces from 2 to 3. At 3, the applied pulse field 
falls to zero. The value 3—3' corresponds to an increased conductivity in relation 


Fig. 6 


to the case of no field. Disconnection of the field at instant 3 corresponds to collecting 
the carriers of opposite sign. As a result, we obtain a rapid drop in volume conductivity 
to 4. The conductivity now slowly increases to 5. 

From the measurements made on polycrystalline layers of PbS, it turns out that 
the slow changes of conductivity corresponding to the segments 2—3 and 4—5 give 
a symmetric wave form. Unsymmetric wave forms are, however, possible, as has been 
observed in the case of CdS monocrystals. 


5. Evaluation of the Time Constants 


The non-steady-state wave forms with which we are dealing in the present measur- 
ing arrangement will not be distorted by this arrangement as long as the numerical 
values of the input circuit time constants are suitably chosen. Evaluation of the time 
constants of the observed wave forms and the time constants of the measuring arrange- 
ment may be represented as follows. 

The rectangular pulse of the electric field applied to the semi-conductor layer 
has a rise time т, = 0.5 usec. We may say in advance that the processes taking place 
in the layer under investigation will not be observed in a time shorter than ту. 

'The time constants of the input circuits of the measuring arrangement may readily 
be evaluated by considering the diagram of Fig. 5 and by assuming that the maximum 
resistance of the cell is 2 MQ. 

Since the output resistances of the measuring bridge are equal to each other, 
and together constitute the matching resistance of the cell, we may use the value 1 MQ 
for our further calculations. The input capacity shunting the grid resistance of the 
first бАК5 valve does not exceed 10 pF. We therefore have here a КС circuit consisting 
of a 1 МО resistance in the bridge arm, a 1 МО grid resistance at the 6AK5 valve, 

and an input capacitance of 10 pF (we ignore the effect of the large coupling capacity 
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10000 pF). The time constant of this circuit is zs== 0.5 x 106 x 10 x 10-1? = 10 usec. 
This value is considerably greater than the pulse rise time, as a result of which the 
leading edge of the pulse is distorted (more sloping) and we cannot observe any change 
in conductivity during this initial time interval. The measurements given below show 
that the time constant of the potential barrier changes are of the order of 1 msec and 
100 msec. Hence the above time constants т, and т; do not influence the shape of the 
observed wave forms. 

Owing to the two time constants examined above, namely the finite rise time and 
the time constants of the input circuits, the wave shape between 1—2 and 3—4 on 
Fig. 6 is not really vertically linear. In addition, from the moment the leading edge 
of the pulse arrives, the occupation of the surface states begins to change which also 
influences the slope of the wave between points 1—2 and 3—4. This is most clearly 
seen by the slight rounding in the vicinity of 2. The true point 2 should be found by 
extrapolating the curve 2—3 to the intersection with a vertical line extended from 
point 1. 

The optimum pulse repetition rate finally adopted for the observations of the 
conductivity changes was relatively low and equalled f = 20 c/s. In this connexion, 
we determined the influence of the time constants of the coupling circuits on the 
pulse amplitude (slope of the top of the pulse). The RC coupling circuit of the input 
stage (Fig.5) consists of a 10000 pF capacitance and a static resistance of 1 MQ. Since 
the stage operates as a cathode follower, its dynamic input resistance is (see Eq. 1) 
Raj, == 20.1 МО = 20 МО. The time constant of the coupling circuit is therefore 
CR = 10000 x 10-1? x 20 x 106 = 200 msec. The pulse width is 1/2. 1/20 sec = 
= 25 msec. Therefore it may be seen that the error in evaluating the rear edge of the 
observed wave form (segments 3—3’ and 5—5”) in Fig. 6) will equal about 12%. 


6. Discussion of the Accuracy of the Method 


The above impulse method was used for microcrystalline layers of PbS for which 
such quantities as concentration, mobility of majority carriers, and minimum density 
of surface states were determined. In order to determine the accuracy of the measure- 
ments, we assume that we had to do with two separate groups of mesurements: 

a. length measurements 

b. electrical measurements 

As regards the second group of measurements, it may be said in advance that the 
instrument errors were not worse than 4%, except for the reading of the values c and d; 
(fig. 7) the latter were determined with an error of about 30%. The greatest error of measu- 
rement was connected with the length measurements, ог more precisely, the measure- 
ment of the thickness Ё of the semiconductor layers. The layer evaporated in the vacuum 
envelope measured approximately 0.1 х 1 x 2 10-5 cm, but could not be measuered 
directly (assuming that the element is not destroyed after the pulse measurement; 
this is important in the case of 'vacuum' layers, since the opening of the envelope basi- 


cally changes its properties). Although the measurements Г = 0.1 cm, m ~ l ет 
can be made with an accuracy of several percent, the measurement of Ё encounters 
great difficulty. А possibility of determining the depth of the microcrystalline layers 
of PbS inside the closed vacuum envelope is provided by interference measurements 
in the visible part of the spectrum or in the near-infrared region. But, however, because 
of the small thickness of the layers and the strong absorption by PbS, this method 
glves only a large approximation. 

There exists another possibility of an indirect determination of the layer thickness, 
namely, by weighing. Unfortunately this requires the opening of the envelope and the 
destruction of the layer. Taking the layer, or its evaporated part on a known area, 
we weigh the layer together with its base and then the base itself. The weight of the 
layer for a known area permits the thickness to be determined. In the case of PbS 
this method gives for k => 1074 cm, an error of the order of 10% (the balance available 
had a sensitivity of 10-5 р). Already for the thickness £ = 2 x 10-5 cm, however, 
this error increased to ~50%. Thus, as may be seen from the above estimations, the 
measurement of the layer thickness decisively influences the resultant accuracy of 
the measurement of the above-mentioned quantities, which in our case, k = 2 x 1055, 
is of the order of 50%. 

It appears that because of the accuracy obtained as well as the necessity of de- 
stroying the layer, it is necessary to work out another method of determining the layer 
thickness. Аз regards the vacuum layers, perhaps it would be possible to work out 
a method of evaporating the layers (with fixed conditions of production: a control 
evaporation) in order to determine the layer thickness from the production conditions 
and the amount of material used for evaporation. 

In the measurements of the density of the surface states, in addition to the difficul- 
ties mentioned already, there is still the question of the potential barrier height (surface 
potential). In our case, this height was not known; only the minimum surface state 
density was determined (this corresponds to the maximunr potential barrier height 
which can be determined). The accuracy of the determination of this value depended 


on the accuracy of the layer thickness measurement. 


Part П 
1. Results > 


The shape of the pulse obtained оп the oscillograph зсгееп is shown in Fig. 7. | 
Ву а simple calculation we obtain the formula for the concentration of the majority 


carriers 


ра KO 1 xs о | 
P 4 (АИ, + ДУ.) kim 1.6 х 10719 Ese! @) 


Because of the capture of the carriers by the surface states, the changes in con- 
ductivity will be smaller; by conductivity modulation we shall mean the ratio of the 
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obtained voltage changes AV, to the total voltage change AV, + AV,, when there 


is still no surface state interaction, viz. 


y АР) 3 
ЛЕНА 9) 


Finally, by the mobility of the majority carriers, we shall understand the quantity 
defined by the formula 


= 20 [ста?/У sec 4 
жены. [У sec] (4) 


Notation 
Е (cm) — layer thickness 
l (cm) — layer length 
m (cm) — layer width 
AV, 
A V, (V) — changes of voltage corresponding to changes of the layer 
AV. conductivity 
AV, 
U_(V) — d.c. battery voltage 
qo (С) — total injected charge 


Preliminary measurements of a series of 24 cells indicated that only some of them 
gave an effect strong enough to permit accurate quantitative measurements. Therefore 


Fig. 7 


four cells, Nos. 8, 10, 21, and 24 were chosen for further investigation. The characteristic 
quantities of the observed pulse were measured, i.e segments a, b, c, and d as regards 
their dependence on the following factors: 


1. pulse repetition rate f 
2. pulse width у 
3. voltage amplitude of the pulse U a. 
4. supply battery voltage U_ 
In addition, also measured was 
5. the decay time т, 
as a function of the pulse repetition rate 


u 
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For each series of measurements the optimum value was determined for the given 
variable which, in the later measurements, played the role of a parameter. The annexed 
photographs of the oscillograph wave forms for the first layer investigated illustrate 
the changes of conductivity. 

The value of the dielectric constant of the glass from which the vacuum envelope 
was made enters into the formulae for the calculations. This value was determined 
from resonance measurements: 


a= HS 


The thickness of the bottom of the glass envelopes was determined from micro- 
scopic measurements. The coefficient of refraction of the glass was determined on 
plane-parallel plates made from the same material. Its value was 


ip == Mr) 


Before making the main measurements, control measurements were made. The 
purpose of these measurements was to see if the connecting of the supply battery U . 
had any influence on the bridge balance. Аз mentioned in the discussion of the method 
of measurement, the bridge with the cell under study is first balanced for rectangular 
electric field pulses. During this time the supply battery U_ is disconnected. After 
balancing the bridge we connect the battery U_ and observe changes in the layer 
conductivity. There exists the possibility that the mere connecting of the battery chan- 
ges the conditions of compensation for the bridge. To check any eventual influence 
of connecting the battery U_, a number of glass envelopes with painted graphite 
electrodes but with no evaporated layer, were examined. Also examined were envelopes 
with shorted graphite electrodes. Finally, layers which did not show any sensitivity 


Fig. 8 


to the applied electric field were also investigated. In all the above three groups of 
control mesurements the bridge balance was not disturbed by the connecting of the 
battery. үй 

Annexed аге several photographs of the typical oscillograph wave forms for 
cell No. 8. . GR 

The results of the measurements for the four cells are given in the table below. 

The results given refer to volume parameters of the semi-conductor layers investi- 
gated. The values of the mobility џ are typical for microcrystalline layers in which 
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the mosaic-structure (internal n-p barriers) lower the mobility of the current carriers. 
The insufficiency of p-type conductivity noted in all four cases is connected with the 
process of producing the PbS layers. During this process, oxygen is introduced and 
becomes and acceptor centre in the PbS lattice (L. Sosnowski 1949). 


Fig. 9 Fig. 10 


Fig. 11 


Table I 


Rl = | -ҰТМТЕІ Ua ana [|] > [2 
= [MOLD cat] ва [mV [шү | —| бее | fose fee] “on [ont i 


5.7 10.15 470x 10121 5.7 x 1074 8.3 
10 | 1.0] 1x10? 28.5 | 5.5 |0.23190х 10-12] 3.8 х 10-4 3.6 
21 | 2.0] 5x10? 18.5 | 14.5 [0.441370 x 1072?! 9.7 x 10-4 31 
24 | 0.2/0.75x 10? 28.5 | 5.5 |0.16580x10-12| 3.8 х 1074 0.5 


2. Discussion of the Assumptions 


In the considerations of the previous sections we tacitly assumed that we are 
dealing with only one type of carrier in the layers — majority carriers, the minority 
carriers being ignored. At present we shall attempt to evaluate the degree of approxi- 
mation involved and improve the results obtained. 

a. Fixing of the values E, and л; 
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If E, is the energy gap at the temperature Т°К, m, and m, are the effective masses 
of the electrons and holes, and m is the mass of the free electron, then the expression 
for the intrinsic concentration of carriers in a semi-conductor has the form 


1 3 з 
па (823 mkT/h?) 2 e-EJ2&T — 4.82 х 1015 х Т? (mere) жама (5) 


Іп order to determine quantitatively the value of the intrinsic concentration it is 
necessary to know two quantities 
Сар expression connected with the effective masses of the holes and electrons 
m 
E, —the energy gap 
L. M. Mackintosh (1956), basing himself on the results of Smith (1954), Avery 
(1954), and Gibson (1952), obtained the following values for PbS: 


Е, =0.3e6eV 
Ао ЖС ЛӨ (уй 


Taking into account, however, the possible errors in determining Ер and the 
ratio b = т,/т, he noted that the value n; is determined to an accuracy of a factor of 5. 
Measurements of the energy gap made thus far (absorption edge measurements) 
gave a value of E, = 0.4 eV. It therefore seems that because of the lack of consistent 
data as regards the value of the energy gap in PbS as well as the possible deviations 
in evaluating п; we may, in our considerations, take the value E, = 0.35 eV. The 
_ value of intrinsic concentration is therefore 1 x 1015 cm 3. 
b. Consideration of the minority carriers 
The semi-conductor conductivity is expressed in the general form 


= neu, + рей» (6) 


where и, and шь are the mobilities of the electrons and holes. The electric field applied 
in a condenser arrangement to the layer acts on both types of carriers, however de- 
pending on the sign of the field, one type will be attracted to the surface and the second, 


repelled. i 
Changes in conductivity with both types of carriers taken into account are expres- 


sed by the formula 


24 
Дс = — P Aneun — E2 Дпець = - = - [nun — Риф] (7) 
| 2 
До = у (n — р) (8) 
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Changes of conductivity Де when the minority carriers are neglected аге therefore 
correspondingly smaller and the value of the majority carriers thus determined is 
smaller than the actual value by approximately the value of the concentration of 
minority carriers which was ignored. 

The actual value will therefore be expressed by the relation (we are dealing with 


holes as the majority carriers) 
р pn (9) 


Actually, p' should be determined from the relation p'n = (p + n)n = n?; from which 
п is found. This is, however, a correction of the second order with respect to n, the 
letter determined from the relation 


п = — (10) 


In connexion with the determination of the corrected concentration of majority сат- 
riers p’ the mobility too undergoes change, viz., 


G G 


йе ^ (pd ye e 


or 


Ес P 12 

Т=П Кт 12) 

The corrections resulting from taking into account the influence of the minority 

carriers are of the order 1% and smaller, and in our case they may be neglected, since 
they lie below the accuracy of the method itself. 


3. Linearity of the Changes of Surface Conductivity 


For the four investigated layers, the volume conductivity and the barrier (surface) 
layer conductivity was found to be of the p type. By introducing hole pulses into the 
investigated layers, an increase in instantaneous and quasi-stationary conductivity 
was obtained. 

There exists the possibility that the potential barriers obtained bend downward, 
and that the type changes to the opposite by the introduction of barriers of the opposite 
sign, in our case n-type. If the number of carriers introduced is increased we should 
obtain, for some value of the field a minimum barrier layer conductivity (minimum 
quasistationary value). 

In his work on p-type germanium Brown (1953) showed that a minimum surface 
conductivity is obtained for sufficiently strong electric fields applied to the surface 
of the layer. The obtaining of a minimum conductivity supplies information about 
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the surface of the semi-conductor, which allows the surface potential (potential barrier 
height) of the sample in a free state (without the electric field) to be determined. 

In connexion with the above possibilities, a layer of PbS specially evaporated on 
mica (in a vacuum envelope) was investigated from this standpoint. This arrangement 
made it possible to apply a considerably stronger field than previously used. The 
investigated layer showe a quasi-stationary p-type conductivity. N-type carriers were 
injected into the layer. 


The results of the measurements are shown in the table below 


Table 11 
ТМ Б? ао 4 TE PES 40 100 
AVe mV 19.8 37.0 74.2 187.7 388 736 2000 
AVc[U | шУ/У 19.8 18.5 18.5 18.7 19.4 18.4 20 
U „n — is the amplitude of the rectangular pulses applied to the layer in a condenser 


arrangement 
Аз it results from the above data, the expected minimum surface conductivity 


does not occur — the changes of conductivity in the limits of the applied fields 0.5— 
50 КУ/ст, are linear. 


4. Determination of the Surface State Density 


The problem of the space charge region was discussed earlier by many authors 
(Schottky, Mott). In our considerations in connexion with determining the wave form 
of the electrostatic potential in the vicinity of the semi-conductor surface, we made 
use of the formulations of Garret and Brattain (1955). The terminology found below 
will be used in the following sense: 
potential barrier (surface potential) — difference between the electrostatic levels 
inside the semi-conductor and just below its surface 
barrier layer (surface layer) — the layer in which the barrier potential (changes of 
electrostatic potential) occurs = 
surface charge — charge in the surface states 
barrier charge — charge in the barrier layer 
Schottky barrier layer — the barrier layer whose conductivity is of a type opposite 
to the volume conductivity ; this layer is depleted of charge (it is also called the inversion 
layer). 


Notation 


n; — intrinsic concentration of holes and electrons in the semi-conductor 
y — electrostatic potential | 

ро — electrostatic potential inside the semi-conductor 

Po — Fermi potential 


£ 
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The variation of the electrostatic potential at the surface can be determined 


from the Poisson equation 


Др = —4ло[{вв, (13) 


In the measurements that were made. it was difficult, because of technical reasons, 
to use fields stronger than 50 kV/cm. At the same time, the lack of a minimum, and 
even. the lack of a deviation from linearity, suggests that the potential barrier is large 
in relation to kT/e. 

It should also be mentioned that the linear character of the changes of conductivity 
were also obtained for the injection of carriers of opposite sign, that is of the p type- 
The character of the changes in the conductivity pulses, as well as the respective time 
constants, did not undergo a qualitative change. This problem will be discussed below. 

In the general case, we have to do with three regions. The part played by each 
in shaping the potential barrier depends on the value of the parameters of the semi- 
conductor. The quantities determining the shape of the barrier are 


В(Фо— vo) 


Д == РОН роп Jipa == ë and у = f(y — vo) 


It follows from the discussion of the variation of the electrostatic potential that 
in the vicinity of the surface of the semi-conductor it is sufficient, in our further 
considerations, to use a parabolic approximation, that is the Schottky approximation. 

The value of the surface potential U = y — wy, as a function of the thickness of 
the barrier layer, is expressed by the known formula following from Schottky's as- 
sumptions: 


€ 


ya 
x= 1, T y = E (polni) * [B(y — фо) = | 


80 
2леро 


% 
| 0% (14) 


Ав а result 


длеро x2 


U= (15) 


2, 

In accordance with our earlier considerations, the electric field applied рег- 

pendicularly to the surface of the semi-conductor introduces a charge into the interior. 

This changes the way the surface states are filled and the charge distribution in the 

barrier layer. As a result, the surface potential U changes by a value dU. Connected 

with this is an additional charge in the surface states 9, = Ме, where Л, is the 
surface state density рег ст? per volt. 


We assume: 


a. that modulation of the semi-conductor conductivity under the influence of the 
electric field is smaller than unity; 
Ь. that the ratio of the barrier thickness to the layer is smaller than unity. 
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The expression for the surface state density then takes the form 


є % 
ат тоз ЕА ior [У-1 ст] (16) 

It is seen from the above relation for the surface state density N ‚ that among the 
unknown quantities is the surface potential U actually existing on the surface of the 
layer. In this connexion, only the minimum value Л, can be determined if the highest 
possible potential barrier is assumed. The maximum value of U,, will be determined 
separately for each cell because of the different values of p (the position of the Fermi 
level depends on this). 


Table III 
и N; 
сш?/У sec У lem 2 
8 5% 0 8.3 9 2 1012 
10 90 52 1018 3.6 5 x 1012 
ӘЛІ 1152 1025 81 16 >< 1012 
24 9 x 1016 0.5 ІЛ x 1014 
Part Ш 


l. Discussion of the Kinetics of the Barrier 


The measurements of the time constants of the conductivity changes in the layer 
were made by the exponential time base method. This gave the values тү == 1 msec and 
т» == 100 msec. Different physical processes connected with the surface of the layer 
are responsible for each of these values. (Volume processes in the layers cannot play 
any role here, since their time constants are of the order of 1 usec.) From what was 
said previously the investigated layers are of the p type. The barrier layers are of the 
same type, but with a higher concentration of carriers. If we introduce n-type carriers 
into the layer, they are drawn to the surface by electric field, thus upsetting the barrier 
layer — surface state equilibrium of charge. There then must occur a change in the 
filling of the surface states. It appears that such a change taking place through the 
potential barrier would be connected with a change in time of the order one msec. 

As regards the large time constant т, == 100 msec, it seems that it should be con- 
nected with other physical processes occurring on the surface. As is known from the 
work of Brattain and Bardeen (1952), the work function depends on the adsorption 
of ions on the surface. The change of the potential barrier caused by the application 
of the electric field produces a change in the work function and therefore disturbs 
the surface equilibrium of the adsorbed ions. Adjustment of the system to a new state 
of equilibrium would be connected with changes in the number of ions adsorbed on 
the surface. These changes take more time. 


£ 
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The previously mentioned symmetry of the conductivity changes in the case 
of the injection of n-type and p-type carriers requires separate examination. From 
purely qualitative considerations, it may be concluded that in the layers of the type 
investigated (p-type volume, p-type surface, barrier bending upward) there should 
occur a considerable difference in the value of the small time constant ту for the injec- 
tion of n or p carriers. In such a case, the introduced electrons must occupy surface 
states by going through the potential barrier, which would extend the time constant 
for the attaining of equilibrium, while in the case of holes, no such obstacle exists. 
Quite the opposite, the potential barrier for the holes, is strongly reduced and the 
changes in the occupation of the surface states in the case of a hole charge in the 
barrier, should occur in a considerably shorter time — of the order of usec. In the 
meantime, no asymmetric processes are observed. In our treatment thus far, we did 
not take into account the influence of the base (mica, glass) on the surface properties 
of the layer. It should be supposed that the potential barrier on the surface 
of the layer is more complex than assumed. It appears that the double-barrier model, 
as shown in Fig. 8, should explain the observed symmetry of the conductivity changes 
(R. H. Kingston 1956). A barrier of this type (E. T. Goodwin 1939) now makes impossi- 
ble a direct exchange of charge between the holes and the surface states (the tunnel 
effect being neglected). 

The only possible way of changing the occupation of the surface states must be 
through the resultant barrier, without regard to whether n-type or p-type barriers 
are introduced. From this it follows that such a barrier would symmetrize the process 
of attaining surface equilibrium. 


2. Applicability of the Pulse Method 


The pulse method worked out for investigating semi-conductors permits many 
parameters to be determined by means of one observation. This is essential where 
successive measurements made by two different methods can change the properties 
of the investigated object. The range of applicability and usefulness of this method 
for investigating monocrystals and polycrystalline materials is different. As regards 
the investigation of monocrystals, it should be said in advance, that there exist many 
methods of measurement connected with transistor technique which allow the determi- 
nation of such quantities as concentration and mobility of current carriers in a simple 
and accurate manner (for monocrystals the accuracy of the dimensional determina- 
tions is greater). On the other hand, it seems that this method creates new possibilities 
for investigating the state of the surface of semi-conductor monocrystals. 

It should be mentioned that, independently of the author, С. С. E. Low used an 
analogous pulse method to investigate germanium and silicon monocrystals (1955). 
As should be expected, the conditions of measurement for monocrystals are much 
easier because of the smaller intrinsic concentration of current carriers and because 
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of the greater current mobility. This permits much greater changes of conductivity 
to be obtained with weaker compensation in the bridge circuit. 

As regards the investigation of the layer, the only generally used method of de- 
termining the concentration and mobility is the method based on the Hall effect. This 
method requires, however, a specially prepared layer. Additional electrodes have to 
be applied to the layer. These are necessary because of the voltage measurements 
connected with the Hall effect. The difficulties in compensating the circuit, as well 
as the small operativeness of this method, leads one to believe that the pulse field 
method will successfully compete with the method based on the Hall effect, as regards 
the determination of the quantities in question. 

The method described in this paper for investigating semi-conductor layers was 
first reported at the Physicists Conference in Spala in September 1954. The results 
of the first measurements made on microcrystalline layers of PbS were published in 
a paper written by Professor L. Sosnowski and the author in the Bulletin de l'Académie 
Polonaise des Sciences, Cl. 3, 2, 101 (1955). 

I would like to express my particular gratitude to Professor L. Sosnowski for 
his choice of the topic, his direction of the work and concern in its progress. 

I thank Dr. T. Piwkowski for preparing the vacuum cells with the lead sulphide 
layers. I also thank Mr. J. Witkowski for making mechanical arrangement used in 
investigating these layers. 


KPATKOE СОДЕРЖАНИЕ 


Г. Ржевуский, Исследование слоев сернистого свинца посредством метода электри- 
ческого импульсного поля. 


Исследован ряд микрокристаллических слоев PbS, напаренных в вакуум- 
ных конвертах. В. настоящей работе шире разобран, раныше разработанный 
(Г. Ржевуский и Л. Сосновский 1955) метод исследования полупроводящих слоев, 
состоящий в применении электрического импульсного поля, прилагаемого к слою 
конденсаторным способом. 

Опираясь на вышеуказанный метод установлены: 

а) тип изучаемого слоя, 

b) тип барьерного слоя, ~ 

с) концентрация большинственных носителей, 

а) их подвижность, 

е) минимальная густота поверхностных состояний, 

f) определены постоянные времени релаксации барьера, 

g) продискумирована кинетика изменений проводимости борьерного слоя. 

В заключение обсуждена точность импульсного метода в применении к изу- 
чаемому роду полупроводящих слоев, а также всякие возможности его приспо- 
соблений. 
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The excitation energies of the ethylene molecule were calculated using in the interac- 
tion integrals of the m-electrons wave functions depending explicitly on the mutual electron 
separation. The resulting energies are in satisfactory agreement with experimental values. 


1. Introduction 


It is well known that one of the main difficulties still to be overcome in the mo- 
lecular orbital theory is to calculate the electron correlation due to the electrostatic 
interaction between electrons. If the wave function of the electrons is represented by 
a determinant built up from one-electron molecular orbitals, then this wave function 
accounts for the correlation of electrons with the same spins. It is the exclusion principle 
which keeps these electrons apart, and the electrostatic interaction plays only a secon- 
dary role in their correlation. However, in a pair of electrons with opposite spins the 
position of one electron relative to the other is influenced only by the electrostatic 
force between them. The wave function for these electrons is symmetrical in their 
space co-ordinates and does not prevent both electrons from being in the same place 
at the same time, which, of course, cannot happen because of their electrostatic repul- 
sion. Therefore, to improve the theory we must use for pairs of electrons with opposite 
spins wave functions which represent the tendency of these electrons to keep apart. 
| The problem of the relation between two electrons with opposite spins and occu- 
pying the same molecular orbital has been discussed by various authors (Coulson and 
Fischer 1949, Mueller and Eyring 1951, Lówdin 1953, Itoh and Yoshizumi 1955, 
Lennard-Jones and Pople 1951, Hurley, Lennard-Jones and Pople 1953), and several 
attempts have been made to generalize the molecular orbital theory in such a way that 
it will take into account the electron correlation. James and Coolidge (1933) showed that 
the electron correlation plays an important role in the stabilization of the molecule. 
They obtained for the binding energy of the hydrogen molecule a value which was in 
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excellent agreement with experiment, and they were also able to show that it was 
possible to obtain this only when using a wave function depending explicitly on the 


mutual electron separation aS = |= rÀ To obtain agreement with experiment, 
however, they had to use an elaborate wave function with many variational parameters, 
and there is no hope at present that this method can be applied to more complex 
molecules. Wave functions depending explicitly on the mutual electron separation 
have been used only for the hydrogen molecule. 

Nevertheless, it seems that it is neglect of the correlation which is responsible 
for the failure of the molecular orbital method of calculating molecular energies, and 
the purpose of the present paper is to suggest a method which to a certain degree 
accounts for the correlation of electrons and which can also be used for calculating 
the excitation energies of conjugated hydrocarbons. We shall apply this method 
first to the ethylene molecule, for which, in spite of its simplicity, the results of the 
conventional molecular orbital theory are not in agreement with experiment. In 


a subsequent paper we shall apply this method to a more complicated case, viz., to the 
benzene molecule. 


2. Тһе Two-Electron Approximation 


In the molecular orbital theory of conjugated hydrocarbons we construct for 2n 
z-electrons of the molecule the one-eletron wave functions фу, d»... Jo, called mole- 
cular orbitals, which are linear combinations of atomic orbitals. In the ground state 
of the molecule, n molecular orbitals corresponding to the lowest energies are doubly 
occupied by electrons with opposite spins. The wave function for the ground state 
of the molecule is then represented in the form of a determinant 


d. (1)«(1) ADLA) al)... | 
p- | AD) $050) &0)«09).. 
У») (904) BBG) HB- 


ACEN 477% 9 s» e» wre a ae 


.-. 


thus satisfying the exclusion principle. 


The wave functions corresponding to the excited states of the moleculeare obtained 
by lifting one or more electrons from the orbitals ф, (k < n) to ће orbitals ó, (12 п). 
If one electron is lifted to a higher level then there are two electrons which are not 
paired. Their spin function may be either symmetrical or antisymmetrical, correspon- 
ding to the triplet or singulet state respectively. In the first case the wave function 
must be antisymmetrical in the space co-ordinates of the two electrons, and thus, 
owing to the exclusion principle, takes into account their correlation. Only the electrons 
with opposite spins and occupying different molecular orbitals, as well as those occupy- 
ing the same molecular orbital, are not correlated in this theory. 
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Let us consider first the singlet states and assume for pairs of electrons with 
opposite spins the two-electron wave functions y,(r,7r) which include the correlation 
of the two electrons. Then the many-electron wave function built up from these two- 
electron functions so as to satisfy the exclusion principle will account for the correlation 
of all electrons in the molecule. This method was first proposed by Fock (1950) and 
Hurley, Lennard-Jones and Pople (1953). Fock has also given a method of building 
up from two-electron and one-electron functions a many-electron wave function corres- 
ponding to electrons with a given resulting spin. 


For the singlet states this wave function may be written in the following form 
Y= (ai) > 92 (—1)” (—1)° РО т (71:72) Po (78, ee V (u 12412 (2) 
ESO 


where #,(7op 41) 72, 49) is a normalized wave function symmetrical in the space co-ordi- 
nates of the two electrons, P is the permutation operator for electrons labelled with odd 
indices, and Q is the permutation operator for electrons labelled with even indices. 


The energy corresponding to the wave function (2) can be calculated from the 
formula 


На (3) 
where 
1 1 
E н dest (4) 
i 47 
апа 
l5. ) j 
Hy =— > Vit ү; (5) 


V; in (5) represents the interaction potential between the i-th electron and the nuclei. 
In the case of the conjugated hydrocarbons V; includes also the interaction potential 
with the o-electrons. 

Assuming for the two-electron functions the orthogonality condition 


f Palis rj) V, (ri Ta) dv; = 0 for n Z m, 7 (6) 


which imposes restrictions on the wave functions y; (r; r)) which may be used, the 
energy of the system may be expressed in the form 


Е = У) [2(p | Alp) + (рр|С|рр)] + > [2(pp| С|чя) —(рд|С\|рд), (1) 
b PES A 


where 


(p |H |р) = ff Pp (т гә) H, Wp (ть) т, Ат, 


| * 1 
(РР | С | рр) = ПЕ (r 72) та Pp (ту, т) dt, ать, (8) 
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; * * I 
(pp | С | qq) = ni (71, ra) varo, ro) T (7, T4) фа (ға, ғә)4ту dz, тз Ата, 
13 


ж ж 1 7 
(pq | G | pg) - [||| = (гі, ra) ра (ra, r2) TTA Фь(ға» ra) Va ("1s го) дт, dt, ат» Ата. 
) 13 


Тһе terms (р|Нір) іп (7) represent the kinetic energy of the electrons and their 
interaction with the core (nuclei and o-electrons), (pp|G|pp) represent the interaction 
between the two electrons in the same pair, (pp|G|qq) and (pq|G|pg) are the Coulomb 
and exchange parts of the interaction energy between two electrons belonging to diffe- 
rent pairs. 

Formula (7) was first derived by Hurley, Lennard-Jones and Pople (1953). When 
V, (rj т) = Bn (ғ) ó, (г), Le. in the one-electron approximation, we get from (7) the 
well known expression for energy obtained by Fock (1930). From (7) equations can 
be derived which correspond to the Fock equations in the one-electron approximation. 

Hurley, Lennard-Jones and Pople (1953) were able to get many interesting results 
from formula (7). However it is not possible to use this formula, without simplifying 
it, for our purposes. The calculation of the integrals representing the interaction of 
electrons belonging to different pairs would be very difficult if the wave functions 
depending explicitly on the mutual electron separation were: used in the calculations. 
It seems, however, that the most important difference in the contributions to the 
total energy given by formula (7) as compared with those in the one-electron approxi- 
mation will appear in the terms representing the interaction of electrons in the same 
pair. This is why we propose to account for the correlation only in the (рр|С|рр) 
terms and to neglect it in the others. A discussion of this assumption will be given 
below. 

Using the two-electron functions for pairs of electrons with opposite spins, we 
may write the wave function for the triplet state in the form (Fock 1950) 


d (Д) ЕЕ) 22 > (--1)”(--1)2 РО ж (ғ), т)... 
Yn—1 (Топ—з› Tən—2) Pr (Ton —1) Фі (Ton) - (9) 


P denotes here the permutation operator for electrons labelled with odd indices in 
the two-electron functions, Q denotes the permutation operator for electrons labelled 
with even indices in the two-electron functions and for the electrons which occupy 
the one-electron orbitals. 

Using (9) we can calculate the energy for the triplet state, and in the approximation 
proposed above, the lowering of the triplet state energy due to the correllation will 
be determined by the (pp|G|pp) integrals which represent the interaction of two 
electrons having opposite spins and belonging to the same pair. In the triplet states 
of a2 n-electron system there are n—1 pairs of electrons with opposite spins, and hence 
п — l integrals of the (pp|G|pp) type in the energy expression. 
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The Ethylene Molecule 


Let us consider now the ethylene molecule, for which a detailed molecular orbital 
investigation was performed by Parr and Crawford (1948). From the two carbon 2px 
atomic orbitals p, and y, we build two molecular orbitals 


fı = (2 + 2 Sere (9, + Ф), 
2 = (2—25) (p, — Фу), (10) 


where 
Se fo: (r;) Ф, (т) dv; 11 
is the usual overlap integral. 


The ground state (N) of the ethylene molecule is then represented by the wave 
function 


Dy = d, (ту) d (ro), (12) 
the excited singlet states (V, Z) Ьу the functions 
Фу = 27% [ф, (r) % (тә) + $, (тә) > (r1). (13) 
Dz = ds (ту) $» (тә), (14) 
апа for the triplet state (Т) we have 
Dr = 2-% [ds (11) ds (т) — d (ғ) Po (r)]- ! (15) 


The energies corresponding to the four wave functions (12) — (15) are given 
by the following expressions 


Ey = 21 + Ju: 
Ez = 2, + Л», (16). 
Ey =1,+ dc Лә t+ К, 
Er = h + h + Je — К» 
where T, Jmn and Kj, are defined by 


I, — |6) G) H. é, (ту) йт, 
Jmn = j [| Dn (r ó, (r 2) = #ы(т)фа(тә) dv, тз, (17) 


Ky = | [| di (r) ф (тә) - фу (rs) б» (r) d, d. 


Only the wave function Фу. includes the correlation of the electrons. In the 
singlet states the electrons move independently of each other in the field of the molecule. 
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То take into account the correlation also in the singlet states we shall use instead of 
Ф, (i = N, Г, Z) the functions 
V, = N, Ф, (т, ту) тә (18) 


where N; is the normalizing factor, which can be calculated from the formula 
es [Г D; (ry. ть) Ф, (ry тә) (ға)! d Gre ae (19) 


The wave function (18) is the simplest one depending explicitly on the mutual electron 
separation. This function vanishes when both electrons are at the same place, similarly 
as in the case of the wave function Ф. for the triplet state. 

Assuming now that we may neglect the effect of correlation on the terms represen- 
ting the kinetic energy of the z-electrons and the energy of their interaction with the 
core, we shall have the following expressions for the energies of the singlet states 

Bi se ly ЕЛ. 
E; = y+ Jp, (20) 
fy=h+h+Jy 


where 
* 1 
12 


: The integrals over molecular orbitals expressed in terms of integrals over atomic 
orbitals will be given by 


(аа | т, | aa) + (aa | rg. | bb) + 2 (ab | ть | ab) + 4(аа | ть» | ab) 
(аа | гу? | aa) + (аа | T12? | bb) + 2 (ab | ть? | ab) + 4(аа | т1о? | аЬ)? 


y wa (аа | г | aa) + (аа | r, | bb) + 2(ab|r, | ab) —4 (aa | rio | ab) 
(aa | тү? | aa) + (аа | r5? | bb) + 2 (ab | rj? | ab) — 4 (aa | rs? | ab)’ 
- (aa | ee | aa) — (ab | rig | ab) 
des (аа | ry” | аа) — (ab | т? | ab) ° с 


where 


Jn = 


(abs lod) = | өп) 92 (г) (ғә) v, (r) Фа (тә) ат, йт, (23) 


Q; denotes a 2pz atomic arbital of the i-th carbon atom. Jy and J, correspond to Лу 
and Ja in the one-electron approximation, and the integral Jy to ЛЬ + Kyo. 

The values of integrals over atomic orbitals of the type (аі (са) and (а са) 
can be found in the tables of Kotani, Ishiguro, Hijikata, Nakamura, and Amemiya 
(1953), and the integrals of the type (ab|r;j|cd) and (ab|r,,2|ed) can be evaluated in 
terms of the auxiliary integrals tabulated partly by Kotani, Amemiya and Simose (1938). 
In the calculations of some (ab|r,,|ed) integrals, however there appear auxiliary integrals 
which have not been tabulated. We have estimated the values of these integrals by 
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using the formula 


(ab|r,əled) = [(ab|ed)(ab|r,,?|cd)]^ ; (24) 


(ab|cd) is, of course, the product of two overlap integrals S,, and 5,4. This formula 
applied to the (ab|r,ə|ab) integral gives the value 0.235, which differs by only 6 per 
cent from the exact value 0.221. The values of all integrals over atomic orbitals are 
given in Table I. They were calculated by assuming the usual value of the effective 
charge Z = 3.18, and the internuclear distance R — 1.33 А. The estimated values 
are marked by the asterisks. 


Table I 


Integrals Over Atomic Orbitals 


n--—1 = 0 nl ncm 
(аа | г," aa) 0.6238 ` . 1.0000 2.192 5.906 
(aa | г12"| bb) 0.3441 1.0000 3.402* 12.204 
(аа | г," ab) 0.13705 0.2869 0.845* 2.488 
(ab | ғ," ab) 0.04214 0.08233 0.221 0.670 


* denotes estimated values. 


Using the values of the integrals given in Table 1, we calculated the energy depres- 
sion resulting from correlation of the electrons: 


App а етеу 
AE, = Ja + К, = 1:00 eV 

AE, = Jon —J7= 5,63 eV 

AE, =0 (25) 


The resulting excitation energies are given in Table II, where they are compared 
with experiment and with the results of the standard molecular orbital method. 


R Table II 
Electronic Energy Levels of Ethylene (eV) 


State MO + configu- | MO + correlation, | MO + correlation Expeinental 
ration interaction!| function (26) function (18) 

N 0.0 

T 6.4? 

V 1.63 

Z 


1 Parr and Crawford 1948. 
? Potts 1953. 
3 Wilkinson and Mulliken 1955. 
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Table Ш 
Inteerals Over Molecular Orbitals (eV) 


MO semi- MO + corre- 
empirical? lation 


Integral MO! 


JN 13.14 8.79 
jp 17.24 10.24 
Jz 13.48 7.85 


1 Parr and Crawford 1948. 
? Pariser and Parr 1955. 


The energies were also calculated by using, instead of V ; = М, Ф, (гг) гә, the 
wave function 


Y; = Ф; (rj 79) + AY (гі, гу), (26) 


where 2 denotes a variational parameter. There was little hope, however, that the 
variational procedure could give any satisfactory results, since we assumed for the 
energy terms corresponding to the kinetic energy and to the energy of interaction of 
the z-electrons with the core the same values which were obtained by using the con- 
ventional molecular orbital method. The excitation energies resulting when the func- 
tion Y; is used in the calculations are also given in Table II. The corresponding values 
of the variational parameter 2 are: — 1.439, — 1.019, — 4.845, for the N, Г and 
Z states respectively. 


4. Discussion 


The agreement with experiment which we obtained using the wave function of 
the form У = N, Ф, (г, го) Гір is even better than could have been expected on account 
of the approximations which we introduced in the calculations. No doubt this agreement 
is partly fortuitous. In the triplet state we considered only the correlation arising 
out of the exclusion principle. The wave function which we used for the singlet states 
is also not the best one to consider the electrostatic. correlation of the electrons. 
А better, but more complicated, wave function would certainly give smaller values 
of the interaction integrals. However, it was shown by Lennard-Jones and Pople (1951) 
that considering the effect of correlation the lowering of the repulsive potential energy 
between electrons is partly offset by a slight increase in the kinetic energy. It may be 
that the errors in the kinetic and potential energies cancel each other. It seems that 
wave functions of the type suggested in this paper may be very successful also in 
other similar calculations. 

We see from the results given in (25) that the lowering of the V state energy due 
to the correlation is much greater than that of the N and Z state energies. The reason: 
for this difference is easy to understand. Moffitt and Scanlan (1953) pointed out 
that the V state of the ethylene molecule is ionic in character, while the N, Z states 
are only partly ionic, and the triplet state is totally atomic. They suggested that the 
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failure of the molecular orbital method of calculating the energies of excited states 
is due to the fact that this method is not able to give a satisfactory value of the difference 
between the energy of two isolated carbon atoms and that of their positive and nega- 
tive ions. 

Mulliken (1953) also discussed the N, T, V, Z states of the hydrogen or ethylene 
molecule. He suggested that the results of the molecular orbital method would be impro- 
ved if the varying degrees of ionic character of the molecular orbital wave functions 
for various states were to be taken into account. We could do this by using, in the terms 
corresponding to the ions, atomic orbitals with effective charge Z smaller than in the 
orbitals for neutral atoms, or, still better, by using self-consistent atomic orbitals for 
the ions. Mulliken expected that in this way the singulet states of the molecule should 
be lowered (the V state about twice as much as the N and Z states), and the triplet 
state, which is atomic in character, should remain unchanged. Following these sugges- 
tions Оһпо and Itoh (1955) were able to improve the agreement between the theoretical 
and experimental results. However, it is evident that the same improvement in the 

‘results would be attained by taking into account the correlation of electrons. It is 
reasonable to expect that the effect of the electron correlation is greater in ionic than 
in covalent states; the greater the degree of the ionic character of the wave function, 
the greater the lowering of the energy due to the correlation. These predictions are 
quantitatively confirmed by our results. 

It may also be interesting to compare the results of our method with those Sue 
in the semi-empirical theory developed by Pariser and Parr (1953, 1955). Pariser and 
Parr were led to the conclusion that the failure of the molecular orbital theory in pre- 
dicting the excitation energies is due to wrong values of the Coulomb repulsion integrals 
which are obtained in this theory. In the r»olecular orbital method, the value of the 
(aa|r,.~1|aa) integral is especially overemphasized. The reason for this overestimation 
is easy to understand. This integral pertains to the mutual repulsion of two electrons | 
attached to a positively charged carbon ion, and it can only arrise from the ionic (C7) 
terms in the wave function. So its value should be improved by using atomic orbitals 
corresponding to a C- ion (Mulliken 1953), or, as we would say, by taking into account 
the correlation of the two electrons. Assuming for the (ааг Цаа) integral a value 
obtained from ionization potential and electron affinity, Pariser and Parr (1955) 
calculated the integrals over molecular orbitals and excitation energies of the ethylene 
molecule, which were in agreement with experiment. А comparison of the values 
obtained by Pariser and Parr for the interaction integrals over molecular orbitals with 
the results of our method is interesting. These values are given in Table III and are 
also compared with the results of the standard molecular orbital theory. We see that 
the values of integrals resulting when allowance is made for the electron correlation 
are not very different from those which were evaluated with the use of a semi-empirical 
method. 

A problem which we have not considered in this paper is the effect of the correla- 

tion on the kinetic energy terms and on those representing the interaction of the z-elec- 


£ 
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trons with the core. In the standard molecular orbital method they are usually evaluated 
using the well known Goeppert-Mayer and Sklar (1938) procedure; the detailed analysis 
performed by Stewart (1956) showed that this method is open to several objections, e.g., 
the calculated excitation energies of the ethylene molecule are very sensitive to variation 
of the effective charge Z in the exponent in the hydrogen-like atomic orbitals for the 
L-shell of carbon atom. Stewart’s results, however are based on the standard molecular 
orbital technique and they do not enable us to draw any conclusions concerning the 
method which includes the electron correlation. It can be shown that the correlation 
changes the effect of the variation of Z on the excitation energies. For example, in the 
molecular orbital method a variation of Z from 2 = 3.18 to Z = 2.80 changes the 


T — V interval by about 1.7 eV, while in our method the corresponding change ` 


amounts to only 0.8 eV. However more precise conclusions concerning the effect 
of correlation on the kinetic energy terms, and on the interaction energy of the z-elec- 
trons with the core require a systematic study. 

The author is greatly indebted to Professor.C. A. Coulson, F.R.S., for helpful 


discussions. 


KPATKOE СОДЕРЖАНИЕ 
В. Колос. Корреляция электронов в молекулах. Г. Молекула этилена. 


Исчислены энергии возбуждения молекулы этилена, применяя в интегралах 
воздействия электронов т волновые функции, explicite зависящие от взаимного 
расстояния электронов. Полученные результаты хорошо согласуются с экспери- 
ментальными величинами. 
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The energy levels of benzene were calculated by an extended molecular orbital 
method. In the interaction integrals, between electrons having opposite spins and belonging 
to the same pair, wave functions depending explicitly on the mutual electron separation 


were used. An assignment of the 1B,, «— 14, transition for the experimental band at 6.2 
eV was suggested. 


1. Introduction 


The method used in the previous paper (Kolos 1957, henceforth designated as 
Part I) to calculate the energy levels of the ethylene molecule will now be applied to 
the benzene molecule. The benzene molecule seems of interest from several points 
of view: (i) it is the typical representative of all aromatic molecules, (ii) the results 
of the standard molecular orbital calculations do not agree with those of the other 
methods in their predictions concerning the symmetry of the transitions corresponding 
to the experimental 6.2 eV band, and (iii) the relative location of the first B,,, Bzw 
and Е,, levels depends only on the interaction of the electrons, and is independent 
of the approximations used in the calculations of the energy terms corresponding to 
the kinetic energy of the z-electrons and to the energy of their interaction with 
the core. 


The best results which can be obtained for the benzene molecule with the use of 
the method of antisymmetrized products of molecular orbitals (ASMO) in the linear: 
combination of atomic orbitals (LCAO) approximation are those of Parr, Craig, and 
Ross (1950). The method used by Parr, Craig, and Ross is an extension of the original 
calculation by Goeppert-Mayer and Sklar (1938). Configuration interaction as well 
as all three-center and four-center integrals were included in this method. | 
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268 W. Kolos 


2. The ASMO Method 


In the molecular orbital method applied to the benzene molecule, 6 mutually 
orthogonal molecular orbitals are formed from 6 carbon 2pz atomic orbitals: 


6 
ó, = (6o) 1 У) exp (2лїЕб)р„ [= 0,41, + 2,3 (1) 
k=1 


where p, is the 2рл atomic orbital of the k-th carbon atom. If the interaction between 
the electrons is neglected, the energy corresponding to the orbital ó, depends only 
on the absolute value |Д. Thus the orbitals % and фз are not degenerate, and the orbitals 
фуу and P,a are doubly degenerate. 
For the one-electron molecular orbitals the following irreducible representations : 
of the Dg, group сап be found 
Po C А » фы СОЁ, (2) 
фз C Bs, , фы C Es. 
As the energies corresponding to the orbitals ø, increase with |Д the orbitals 
Фо фу and G_, are doubly occupied in the ground state of the molecule, and the wave 
function for the ground state will be given by 


Фо = фе(ту)фе(гә)фу(гз)фу(т4)ф- ү(г)ф- (т). (3) 


То take into account the exlusion principle ме shall multiply this wave function 
by the corresponding spin function and antisymmetrize it. 

The first excited level arises by lifting one electron from the orbital ó, or f, 
to the orbital ó, or $ ,. Four wave functions belong to this level 


фу = PTP oA 64 (гз) (7а) (ть) (гє), 

Po = фо(ту)фо(тә)ф-1(т)ф-1(т)фу(гь)ф-»(гв)› (4) 
фз = Фо(г1)6о(ғ%)61(ғз)Фі(ға)Ф-а(г5)Ф- (б), 

фа = dor) o(ra)- 1673) (7) (6, (rs) (ro). 


They transform like the products Á, 64, i.e., 


Ф V» Ya Фа C Erg X Ej = В, + Bou + Ey. = (5) 


The linear combinations of the functions (4) belonging to the states of the symme- 
tries B,,, By, and E,, are : 


9 (Biu) = (gi + Ф), 


Вей! 
9 (Bau) = VE (Ф — P2)» (6) 


Ф(Е, 1) = фу, 
P(E) = Фа. 
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When the wave functions (6) are used in the calculations including the electronic 
interaction, the four-fold degenerate level splits up into two single states belonging 
to the B,, and В», representations and one doubly degenerate belonging to the E,, 
representation. 

In a similar way, excited states corresponding to higher energies can be considered, 
e.g., those arising by lifting one electron from the orbital Á, to the orbital 6,5, or 
from фуу to фз, etc., as well as the states arising by lifting two electrons to the higher 
orbitals. However, in addition to the states В;„, В, and E,, discussed above we will 
consider only the state arising out of the lifting of one electron from the orbital фу to 
the orbital d , ,. because for this transition the experimental band at 6.2 eV was assigned 
by Craig (1950) and Parr, Craig and Ross (1950). The wave functions for this doubly 


degenerate state will be 


Ps = Pil) Pr (72) Ф- 13) 1(га)Ф o(rs)ds(re). (7) 
9s = ñ(r1)ó)(r;)ó_1(rs)@_ (га) о(15)6- (т). 


They transform like the product фф, ә, i.e., 


Qs; Ф СА, X Eo, = Ез, (8) 

Therefore 
Ф( Е.) = Фо Ф(Е,,) = Pe: (9) 
The spins of the two electrons occupying the same molecular orbital must be 


opposite, and the corresponding spin function will be «(1)5(2). For the electrons occu- 
pying different molecular orbitals the spin function may be antisymmetric 


a [a(1) 20) — «(2) £N), (10) 
or symmetric 


T [a(1) (0) + a(2) ВО p 


corresponding to the singlet or triplet state, respectively. 
Thus, for the ground state of the benzene molecule we will have the spin function 
of the form 


Хо = %(1)8(2)о(3)8(4)9(5)8(6), | (12) 


and for the excited singlet states 


» = чу OPORBA 1900) — «GL 08 


£ 
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Taking into-account the exclusion principle we will antisymmetrise the product of 
the co-ordinate and the corresponding spin function 


ФО) = ai Уу (—1)РРФ(ГУ)у:. (14) 


where P is the usual permutation operator for the electrons. Then from the formula 
EL) = fou) Н Ф(Гуат (15) 


we can calculate the energy corresponding to the wave function Ф (Г), i.e., the energy 
of the state which has the symmetry Г. The Hamiltonian operator in (15) is given by 


He 2 e (16) 


Tij 


where H; is the Hamiltonian operator for the i-th electron moving in the field of the 
(C H) *% core. 
In this way we obtain 

Е(14,) = 26g + 4; + Yoo ak 8yo + буз — 409 — 264-1, 

E(B,,) = Е, + 6-42 — д — 6 + 20, (17) 

E(B) = E, + 6-12 — дь + £ — 2n, 

EC E,,) = Е, — 0-33 + Sp, 
E(*E,,) = eo + 4&4 + & + 4y gi + Y os E 6y31 + 4y4143— 2091 — 204.,— Oy. —O_12 + до», 
where 
E, = 280 + Зе + & + Yoo + буо + 27 + Зуп + Зу» — 3091 — доз — 44-1 (18) 
£j is defined by 


£j — ТЕДІ (r,)dt, (19) 
and the integrals y;,, Ôj, and 7) аге special cases of a more general integral 
+ 1 
бы? = [| ó, (r) ФА) та Фк Ф") d, ат», (20) 
ово y И (21) 


š = бы", 9 = ба 


The lowest triplet state of the benzene molecule is the ЗВ; , state. The wave function 
for this state can be obtained by antisymmetrizing the product of Ф(В,,) and the 
spin function l 


Дет x^ x(1)B(2)a(3)B(4)La(5)8(6) + «(6)8(5)]. _ (22) 


Ad 
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The corresponding energy of this state will be given by 
ECB) = Е — д 1—65 — 6. (23) 


The integrals over molecular orbitals defined by (19) and (20) can be expressed 
in terms of the integrals over atomic orbitals which were calculated by Parr, Craig, 
and Ross (1950). Taking into accounting the configuration interaction Parr, Craig, 
and Ross obtained a lowering of the energies as compared with those given by the for- 
mulas (17) and (23), and the order of the 1В, and lE, states was changed by this 
interaction. 


3. The Effect of Correlation 


When allowance is made for configuration interaction the molecular orbital 
method includes, to a certain degree, the electron correlation. The ground state energy 
of the benzene molecule was found to be reduced by 2.7 eV by this interaction (Parr, 
Craig, and Ross 1950, quoted by Itoh and Yoshizumi 1955). 

For the ground state of the benzene molecule the effect of correlation was calcula- 
ted recently also by Itoh and Yoshizumi (1955) who used the alternant orbital method. 
For a pair of electrons (for which in the scheme outlined above the same molecular 
orbitals were used) they constructed two different orbitals having large and small 
amplitudes at the alternate atoms. This lowered the ground state energy of benzene 
by 2.35 eV. However Itoh and Yoshizumi did not extend their calculations to include 
also the excited states of the benzene molecule. 

We shall take into account the electron correlation in benzene according to the 
suggestions given in Part I. We suggested there the use, for pairs of electrons with 
opposite spins, two-electron wave functions of the type 


Vn = N, Pn (ri rj) rjj (24) 


where ¢,(7;, r;) is the conventional molecular orbital wave function for the two electrons, 
symmetrical in their space co-ordinates. We suggested also that the electron correlation 


. be taken into account only in these energy terms which represent the interaction of 


two electrons belonging to the same pair. 

Thus we shall follow the method first used by Goeppert-Mayer and Sklar (1938) 
and extended by Parr, Craig, and Ross (1950) which was outlined above, using instead 
of one-electron molecular orbitals the two-electron wave functions. 


Let us denote 
mm = Apr Pn (rj) Pm (r;)/ i) 
Vm. = Е „ОВА + ВДА) for т (25) 


From the properties of the molecular orbitals ¢, it follows that 


ok = V-n-n* L (26) 
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The wave function for the ground state of the benzene molecule built up from 
the two-electron functions (25) according to the method of Fock (1950) will be 


< Pel) = 2 2 у, (CHP РО, alte Vals) (27) 


== | Poo (ro r3) Vii (Ta T) 0-1 in r9) |, 


and for Ше excited singlet states we shall have 
1 
si) = үз {Vool 7) Via (ras 74)0-(75» гв) |+ | Poo (r %)9-і-1(7» 74) P1-2("s> r9) |)» 


il 
PB) = ү? {| Фоо(гі» r2) Vii (ras Га) P-12(75> 76) | ер | Pool" 72) 1-1 (7, T4) Py-2(Ts07@)|}> 


ШЕ = |Poo (т, 72) Vai (73, 74) V-1-2 (or re); (28) 
2 СЕД) = [vii (гі, 72) Ші (ras Ta) Фо» (Te 76)| 

The wave function for the triplet state ЗВ; can be obtained from that for the 1В;, 
state by replacing y.,4 and у, by фф, and $49. , respectively, and then by anti- 
symmetrizing the wave function according to the formula (9.1) from Part І. 

The energy corresponding to the functions (28) can be calculated from formula (15) 
by means of the wave functions Я (Г;) instead of (I7). Following the suggestions 


given in Part I we shall put f (г) = ry in the energy terms representing the interaction 
of two electrons beloging to the same pair, ie., in the integrals of the tvpe 


" 1 
Gmn = JJ mn (ri, rj) s Фит(Г;, rj) dt; ат;, (29) 
ij 
and f(7;;) = 1 in the other terms. We thus obtain 


W (Arg) == Е(14,,) = оо — 2yn + Goo + 263, Ë 
W (hi) = ECR) — yoo — Yu mapasa 9-12 d- Gog бы + 6-і» 
W(B,) = E(B,) — Yoo —Yn — Уз — 0-12 + Goo + Gn + Can 
WQE,) = E(E,,) — Yoo — Yu — yu — дә + Goo + Gn + б» (30) 
W (Ezp) == Е(Е,.) — 2уц — Yor — до + 263 + Со» 
© ЖБ = EQB,) — yoo— Yn + Goo + Gu, 
where Е(Г) are the energies given by formulas (17) and (23), i.e., resulting in thes 


ASMO method before allowance is made for the configuration interaction. The inte- 
grals Gmn were defined by (29). From (26) it follows that 


НИЕ dam. ату 
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Denoting 
м? = [4 (т) ф; (т) тә Pr (тә) ф, (ғә) dt, ат», 
NH = f S 6; (т) é; (r) (ru)° é, (ғ) d; (r9 ат, dz, (32) 
we have 
doe at и ME CO (33) 
and for n Z т 
Ca == (И m—m 4e M, A eese) (NE "ы; == Nous may. (34) 


The integrals over molecular orbitals y;; and д,; which appear in (30) were evalua- 
ted in terms of the integrals over atomic orbitals by Parr, Craig, and Ross (1950). 
The integrals М,%, M, and hence the integrals (С, can also be evaluated in terms 
of atomic integrals. 

If we write for brevity 


єр (a) р a p (0) ра (в) d, d 
== (ab |r” | ed) = (1,1 + u|rg"|l + v, 1 + w), (35) 
where p, is the 2рл atomic orbital of the k-th carbon atom, then by using the formulas 
derived by Parr, Craig, and Ross we can express the integrals М, ы? and №, as linear 
combinations of the 30 different integrals (ab|r45"|cd) 
Мы 6С (ll и Тіҙ |l + v, 1+ w), 
Мы SOC СВ a prgtil + a, 1+ >), (36) 


where 
Cow = (5;[85„) (2 cos [(2zJ6) (iu + kv + lw)] 
+ 2 cos [(27/6) (iu + lv + kw)] + 2 cos [(2zJ6) (ju +- kv + lw)] 
+ 2 cos (2л/6) (ju + lv + kw)] + 2 cos [(2zJ6) (ku +- iv + jw)) (37) 
+ 2 cos [(27/6) (ku + ju + iw)] + 2 cos [(2л/6) (lu + iv + jw)] 
+ 2 cos [(2л/6) (lu + jv + iw)]}, 


in which the "molecular symmetry number" s,, is 1 for integrals with no axis of 
symmetry [e.g., (12 | ri” | 35)], 2 for integrals with an axis of symmetry in the plane 
of the molecule | e.g., (11 | гә” | 14)], and 4 for integrals with an axis of symmetry 
perpendicular to the plane of the molecule [e.g., (11 | ә" | 44)], and the "integral 
symmetry number" s; is the number of distinct members in the set [e.g., s; — 8, 4, 
2 for (12 | гы” | 35), (11 | r^ | 14) and (11 | ғ," | 44) respectively]. 

Some comment is required about the method of calculating the integrals over 
atomic orbitals which appear in (36). The two-centre integrals of the type (ab | rj5?|cd) 
were calculated in terms of the auxiliary integrals tabulated by Kotani, Amemiya, 
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and Simose (1938). Three-centre and four-centre integrals of the same type were 
estimated by the method of Sklar (1939), 1.е., from the formula 


(ab | г," | ed) == 5,9,1 (kk | ry” | ib), (38) 


where S,, and $,, are the overlap integrals, and the centres k and / refer to the hy- 
pothetical nuclei located at the place where the functions p, and рь, ог p, and p, respecti- 
vely, overlap the most. This method, originally used to calculate the usual integrals 
appearing in the molecular orbital method, has been extended by us to include also 
the integrals with (r9)? in the integrand. Two-centre, three-centre and four-centre 


integrals of the type (ab | гр | cd) we estimated by the method proposed in Part I, 
i.e., from the formula 


(ab | ry | са) == [(ab | са) (ab | r12? | са), (39) 


where (ab | cd) is the product of two overlap integrals 5,, 5,4. 

The values of all 30 integrals over atomic orbitals are set out in Table I. They 
were obtained for the effective charge in the Slater-type orbital Z — 3.18 and for the 
interatomic distance R = 1.39 À. The integrals G„„ were calculated with the use of, 
these values and formulas (33) and (34). The results are shown in Table II, where the 
values of the corresponding integrals obtained in the molecular orbital method and 
in the semi-empirical Pariser and Parr's (1953) theory are also given (Си correspond 
to ут and G,, to y, o): 


Table I 


Integrals Over Atomic Orbitals 


ШЕЛ ЕТІП БЕНЕРНЕН Е 1 


(ии) 5.93 2.197 (12 | j^ | 34) 1.461 0.316 
(11; 122) 12.83 3.582 (12 | rj? | 45) 1.814 0.352 
(117133) 26.63 5.160 (12 | r7 | 13) 0.079 0.029 
(11 | 7;;" | 44) 33.54 5.791 (12 | r7 | 24) 0.133 0.037 
(11 | 7;;" | 12) 2.314 0.777 (12 | rj; | 35) 0.187 0.044. 
(11 | rn^ | 23) 4.698 | 1.107 (12 | ri^ | 14) -0.053 0.016 
Q1 | r^ | 34) 7.399 1.390 (12 | rjj | 36) 0.053 0.016 
(11 | r^ | 13) 0.548 0.147 (13 |; | 13) 0.018 0.0053 
(11 ri | 24) 0.715 0.168 (13 | rij | 24) 0.012 0.0044. 
(11 | r^ | 35) 0.852 0.184. (13 | t;;” | 35) 0.017 0.0052 
(11 |n | 26) 0.304 0.110 (13 | r;,” | 46) 0.020 0.0056 
(11 | rj" 1 14) 0.293 0.073 (13 | r^ | 14) 0.006 0.0020 
(11 | n^ | 25) 0.234 0.065 (13 | rj | 25) 0.006 0.0020 
(12 | rj" | 12) 0.570 0.197 (14 | 73" | 14) 0.004 | 0.0012 


(12 | r^ | 23) 0.756 0.227 (14 | rj?" | 25) 0.002 0.0007 
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Table II 
Integrals Over Molecular Orbitals (eV) 


Е MO + corre- MO semi- 
lation empirical? 
6.166 6.82 
6.82 
7.45 
8.06 
7.26 
1 Pariser and Parr (1953). 
2 Parr, Craig, and Ross (1950). 
Table III 


Electronic Energy Levels of Benzene (eV) 


Experimental 
interaction! 
Та 00 0.0 0.0 
se 4.9 
ns 6.2 
S 7.0 
Ls 52 
x í 3.8 


1 Parr, Craig, and Ross (1950) 


The resulting depressions of the energies due to the electron correlation for the 
five singlet states amount to 


E(?A,,) — F id) = C+ yu — Сы = G -- 2.57 eV, 
E(B,) — W(B,,) = C + у + 9-12 — Gaz = C + 3.95 eV, 
E(B,)— W(B,) = C+ yag + 0.,5— С-р = C + 3.95 eV, _ (40) 
Е(1Е,) — W (CE) = C+ уг до Сә = C + 4.06 eV, 
Е Cis) Суво nt e 
+ ógg + Соо — Gn — Со = C+ 4.24 eV, 
where 
C = Yoo + Yu — Соо — Gi (41) 
and for the triplet state | 
_ EGB,) — WGB,) = C. (42) 


The resulting excitation energies are given in Table Ш and compared with 
experimental values and with the results of the molecular orbital method. 
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4. Discussion 


The best results which may be obtained for the benzene molecule by the standard 
molecular orbital method including the configuration interaction are those of Parr, 
Craig, and Ross (1950). The configuration interaction changed the order of the 1B,, 
and 1£,, states and give for the 1B}, state a higher value than for the 1Е,, state. The 
experimental band at 6.2 eV was therefore assigned by Parr, Craig, and Ross for 
a 1E, < 1A,, transition. The first suggestion of this assignment, which is not in 
agreement with the previous results of the molecular orbital method, is due to Craig 
(1950c). Craig calculated the energy levels of benzene using the valence bond method 
and including the polar structures; he obtained for the LE, < 14,, transition the 
value 6.6 eV. This method was a semi-empirical one and the values of the empirical 
parameters were determined (Craig 1950a, b) from the experimental results for ethy- 
lene. However, a very accurate analysis of the ethylene spectrum performed recently 
by Wilkinson and Mulliken (1955) showed that Craig's interpretation of this spectrum 
was rather not the correct one. Therefore the values of the empirical parameters in 
the valence bond calculations should be changed, and this will certainly give other 
excitation energies for the benzene molecule. | 

The influence of configuration interaction on the energy levels is also not free 
from several objections (compare: Stewart 1956) and is certainly overemphasized in the 
molecular orbital theory. For example, in ethylene the non-diagonal element of the 
configuration interaction matrix should be equal to half the difference between the 
first excited singlet and triplet level. The calculated value of this element, however 
is about 7 times greater than would follow from the experimental V — T interval. 
Pariser and Parr (1953), and Stewart (1956) showed that by assuming an experimental 
value for this element the influence of the configuration interaction on the location 
of the ground level becomes very small. 

The excited states of the benzene molecule have been studied also by other authors. 
Міга (1952) considered the influence of the o—z interaction on the energy levels, 
Parriser and Parr (1953) calculated the excitation energies using a semi-empirical 
method, and it follows from the results of these calculations that the experimental 
band at 6.2 eV should be assigned rather for the !B,, < 141, transition and not for 
the 1E,, < 14,, transition as suggested by Craig (1950c) and Parr, Craig, and Ross 
(1950). This is in agreement with our results and therefore this interpretation of the 
benzene spectrum is given in Table III. 

The method used in this paper includes only the electrostatic correlation of 
electrons with opposite spins and belonging to the same pair, and neglects the electro- 
static correlation of electrons with the same spins and of electrons with opposite spins 
belonging to different pairs. The influence of correlation on the kinetic energy of the 
7z-electrons and on the energy of their interaction with the (CgH,)** core is also í 
neglected in this method. 


How does this affect our results? Let us suppose that we use formula (7,I) from 
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Part I, which was derived with the restriction (6,1), and we do not introduce any other 
approximations. Then the kinetic energy and interaction with core terms will have 
the same values for the 1B,,, 1В,, 1Е,, states. Therefore the relative location of these 
states is not affected by the approximations used in the evaluation of the above mentio- 
ned integrals. It was however shown by Lennard-Jones and Pople (1951) that the values 
of the kinetic energy integrals increase if the wave function accounting for correlation 
is used in the calculation. We suggested the use of two-electron wave functions for 
all pairs of electrons in the ground state (and in other singlet states), and one-electron 
orbitals for the two electrons with the same spins in the triplet state. Therefore if 
allowance were made for the correlation in the kinetic energy integrals we should 
obtain a smaller difference between the corresponding terms for the ground and triplet 
state than that obtained in this paper, and our result for the triplet state would be 
improved. 

Let us consider now the other approximations which we introduced in our method. 
The influence of correlation on the energy of interaction between electrons belonging 
to different pairs was neglected. However, the part of energy corresponding to Ё, in (17) 
would have also the same value for the 1B,,, 1B,, and 1Е,, states if two-electron wave 
functions were used in the calculation. The relative location of these states would be 
affected by new values of the exchange integrals corresponding to д, and Ô- which 
appear with the minus sign in the energy formulas (17), i.e., these which were denoted 
Бу (pq | G | pq) in (7,1). Those which appear with the plus sign were included in 
the G „n integrals in our method. However, the relative location of the 1B,, and 1В,, sta- 
Les depends only on the values of the £ and ?) integrals which represent the interaction 
of electrons belonging to different pairs. There are reasons to suppose that the value 
of these integrals would be greater if allowance were made for electron correlation. 
However they appear with opposite signs and so this effect of correlation will be partially 
cancelled. On the other hand, we have £ with a minus sign in the energy formula for 
the triplet state ?B,,. Therefore we may expect that if we were to take into account 
the correlation of electrons belonging to different pairs we would obtain a lower value 
of energy for the 3B,, state. This is the second effect which can improve our result 
for the triplet state.  . yx 2 

For the reasons given above it seems that 18,, — 1B,, interval would not be con- 
siderably affected if we were to take into account the correlation in all energy terms. 
This supports our assignment of the 6.2 eV band for the 1B,,, < 14,, transition. It seems, 
however, that if we were to take into account the correlation in all energy terms, we 
would get for the triplet state ЗВ; а lower value, which is in agreement with experiment. 


KPATKOE СОДЕРЖАНИЕ 
B. Колос, Корреляция электронов в молекулах II. Молекула дензола. 


Исчислены энергетические уровни молекулы бензола, пользуясь расширен- 
ным методом молекулярных орбиталей. В интегралах, представляющих воздей- 
ствие электронов д со спинами противоположными и принадлежащими к этой же 
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о N E 


самой паре, применены волновые функции, зависящие explicite от взаимного 
расстояния электронов. Полученные результаты для сингулетных состояний XO- 
рошо согласуются C экспериментальными величинами. Опытной полосе при 6.2e V 
приписан переход ІВ,, А 
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A model of hot particles in isospace is investigated enabling a physical interpretation 
of charge, isobaric spin, attribute (strangeness), and mass in terms of rotations in isospace. 
This model possesses some analogies with the hydrogen atom: the particle possesses 
a spinning centre, a spinning outer point, an orbital angular momentum, and a size. 
The particle (antiparticle) possesses also a property of a right (left) screw. The spin of the 
centre is interpreted as isobaric spin. Strangeness is replaced by a more physical attribute 
being an extremum component of the angular momentum of the outer point. This angular 
momentum, together with the radius of the particle, determines the rest mass. A privileged 
“third axis“ is no more a property of isospace but of the particle itself. The algebraic 
additivity of the attributes is secured by means of a generalized postulate of contact inter- 
action. Conservation of the total angular isomomentum holds for strong as well as for weak 
interactions between hot particles. А new “particle — isotrino is suggested. 


1. Introduction 


The present situation in the theory of elementary particles is paradoxical as: 
regards the notion of isobaric spint. A great success of this conception in the domain 
of strong interactions suggests the isospin to be not only a formal though convenient 
shorthand but to possess a deeper physical significance. On the other hand, the existence 
of a privileged “third axis“ for electromagnetic interactions and the lack of isospin 
conservation for weak interactions suggest just the opposite. Thus, we are faced with 
the following alternative: If the conception of isospin is to be taken quite seriously, 
we have to expect the existence of a strict conservation law (like the conservation of 
momentum, energy, or charge). Then we have to explain why weak interactions violate 
(or seem te violate) this conservation law, or stressing the lack of isomomentum 
conservation in weak interactions, we have to regard its conservation in the case of 
strong interactions as more or less accidental. Since this last alternative does not 


1 [sobaric (or isotopic) spin will be called henceforth isospin for brevity. We shall also use the names 
isomomentum and isospace. 
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explain anything, we prefer to consider the conception of isomomentum seriously, and 
ask: is the isospin identical with the total isomomentum of the system? Ií it is not, 
we may save the law of conservation of theisomomentum in spite of a lack of conser- 
vation of a part of it — the isospin. To this end we have to generalize the present 
formalism by introducing (besides the isospin) at least one more quantity of a similar 
transformation character. Moreover, we have to reinterprete the privileged direction 
(third axis) to be no more an intrinsic property of isospace but rather a manifestation 
of an anisotropy of the particle itself. 


2. Redefinition of strangeness (Attribute) 


Let us start with a discussion of the well-known formula of Nishijima, Gell-Mann 
and Pais 


b 
о 


where 4 is the charge (in units of the elementary charge), tą is а component of the 
isospin, s is the “strangeness“, and b is the number of baryons (heavy fermions) minus 
number of antibaryons. This phenomenological formula is not unique. The classifi- 
cation of interactions into weak and strong is equally well possible if we write it more 
generally 


q = + із + % + Bb (17) 


where а, В are arbitrary numbers and the + sign accounts for our freedom of choosing 
a positive direction of the “third axis“. A fortunate choice of the factors may be helpful 
for guessing a deeper physical meaning of the connection between charge, isospin, 
and baryon number, whereas a bad choice may abscure it completely?. Deciding 
for a minus sign of tg, a promising choice seems to be В = —a = 1. The “strangeness“ 
defined with this choice of factors will be denoted henceforth by a and called attribute. 
Thus, we have a new formula, wthout the „опе halves“: 


b—q-— tg +a. (2) 


With this choice the attribute assumes the values given in. Table 1. 
The rule of Gell-Mann holds in the following form: for 4a = 0 interactions are strong, 


is Р 
{ог |Да| = зу interactions are weak, for |Да| > 1 interactions are very weak ог none. 


The columns in Table 1 are arranged so that the neutral pion constitutes an axis 
of symmetry. The other particles situated symmetrically with respect to this axis are 


1 
2 Some authors put « = m В = 0 whereby a symmetry is brought about between the baryons № 


and 5 (their *strangenesses" become 1 and —1 respectively). However, it seems doubtful whether such 
a “symmetry” has anything to do with reality on account of a big mass difference between these particles. 
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Table I 


charge conjugate one to the other. We notice the following regularities of the present 
scheme: ! 

(a) Calling all particles to the right of the neutral pion in Table 1 “particles“ and 
those to the left “antiparticles“, we notice that the attributes of all particles are positive, 
those of all antiparticles negative. 

(b) The pion assumes a privileged position since it is the only particle with a vanish- 
ing attribute. 

(c) The attributes of all particles with integer isospin are integer, those with half 
integer isospin are half integer. 

(d) The absolute value of the attribute increases with the mass of the particle’. 

The above mentioned regularities seem remarkable. The fact that the sign of the 
attribute is connected with the existence of two classes: that of particles and that of 
antiparticles, suggests that the physical meaning of the attribute is closely related to 
the particle-antiparticle character. In our scheme the pion assumes (by its a = 0) 
an exceptional position. This seems very satisfactory since the (neutral or charged) 
pion is actually distinguished by the fact that it is the only hot particle for which the 
operations of charge conjugation and charge symmetry coincide. Moreover, the neutral 
pion is the only particle identical with its antiparticle. Finally, the pion is distin- 
guished among all hot particles by an exceptionally small mass. 

The remaining regularities inherent in our scheme will be discussed in more 
detail in the next sections. 


3. Interpretation of attribute as a third component of an isomomentum 


The fact that all particles with an integer (half integer) isospin possess integer (half 
integer) attribute seems also remarkable and suggests an intimate relation between attri- 
bute and isospin. The nature of this relation may be clarified by taking into account that 


з For the traditional choices of the parameters х, В these regularities were lacking. E. g., assuming 
В = 0, the signs of a for N and Z are opposite notwithstanding the fact that neither is an antiparticle. 
Thus, any relation between the sign of strangeness and the particle-antiparticle character was lacking. 


Е 1 ` L ` 
Moreover, by putting В = 0 [ or В = i) not only the pion but also A (or N) acquired a vanishing 


strangeness. This suggested an analogy between л and A (or N). Obviously no such analogy exists. 
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the sum ts + a isa strictly conserved quantity for strong as well as for weak interactions 
(involving only hot particles listed in Table 1). This fact may be interpreted in the 
following way: not only tg but also a denotes a “third component“ of an angular isomo- 


1 : 
mentum іп (the same) isospace. А change of t, by At, = > does not necessarily mean 
a violation of isomomentum conservation but may be interpreted as an exchange 


: IY. 1 $ 
between two constituent parts of the system: t being diminished by >: athird compo- 


nent of another isomomentum is increased by the same amount Ла = + , or vice versa. 


2 
Thus, the law of conservation of the third component of isomomentum may be secured 
by assuming that (the hot) particles possess a structure in an isospace, i.e. are composed 
of two parts, one endowed with an isospin, the other with another angular isomomentum 
whose third component is identical with the attribute a. 

Thus, the hot particle structure (in isospace) is analogous to the atomic, structure 
(in the usual space), its total angular momentum being composed also of two parts: 
the spin of the nucleus and the angular momentum of the electronic shell. 

The fact that both t and а are simultaneously integer or half integer may be 
stated in the following way: hot particles are always Gas , l.e. possess an integer 
total angular isomomentum. 

The law of conservation of charge (or, strictly TE of the difference b—4) 
may be traced back to the conservation of (the third component of) the total angular 
isomomentum, and the quantity 5—9 may be interpreted as a rotation in isospace. 


4. Connection between mass and attribute 


Table 1 shows that the mass increases with the absolute value of the attribute 
and assumes an exceptionally small value for а = 0. This fact suggests a relation between 
mass and absolute value of the attribute. 

А widely spread opinion is that the rest masses of elementary particles are expli- 
cable in terms of a self-action and therefore cannot be predicted unless a consistent 
theory of interactions will be formulated. We are not convinced by such an argument. 
It seems that an explication of such big masses as those of baryons in terms of field 
masses would require very large coupling constants. This seems to contradict recent 
results in meson theory where the nuclear forces coupling constant appears to be 
considerably less than unity (g? — 0.08). With such g? the ratio 9/265 of the electro- 
magnetic to the remaining mass of the pion is still of the order e?/g? so that the whole 
mass of the pion is explicable as a field mass. But this ratio does not agree with the 
corresponding mass ratios in the case of baryons. Their rest masses seem to be too 
large to be a mere results of a self-action due to nuclear forces. Therefore it seems 
reasonable to look for quite a new explanation of the heavy particle masses regardless 
of whether a satisfactory theory of nuclear forces is available or not. 


Screw Model of Particles : 283 
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Let us assume, as a working hypothesis, that the main part of the heavy particle 
masses is not a field mass (in the traditional meaning of this word) and call it “mecha- 
nical mass“. This mechanical mass is zero only in the case of the pion (leptons are 
disregarded). It seems remarkable that the mass ratios of the baryons are very close to 


ттд mp 0 йс = l g RE 28 0а (3) 


The probability of an accidental coincidence of these numbers. with the experimental 
mass ratios is less than 10-3, so that is seems reasonable to expect a theoretical expla- 
nation of the “magic regularity as suggested by (3). 

According to the individual values of the attribute a given in Table 1, the ratios (3) 
may be written in the form 


1 \% 11% 
MN e mA Та: ТЕ = ics (ala + Dyk (ala D(a + 1) . (4) 


(4) not only reveals that the rest masses are connected with (the absolute value of) 
the respective attributes but also that they are (simple) functions of the length |a| of 
an angular momentum vector. It is striking that the masses of A, 27 (i.e., of particles 
with integer isospins) are connected directly with the square of an angular momentum 
(with the eigenvalues a(a + 1), where a = 0, 1,... is the length of the angular mo- 
mentum vector), whereas the masses of N, & (1.е., of particles with half integer isospins) 


: : А . : 1 
are connected directly with the operator j of Dirac (with the eigenvalues а + — , where 


2 
а = =, = ,... is the length of the angular momentum vector). 


A similar connection seems to hold also for mesons: 


Ча 
my : т, = (а (а -1))^: (« + A 5 (5) 


which accounts for the fact that the mechanical mass of the ріоп vanishes (а = 0) 
whereas that of the meson K is different from zero,-in a qualitative agreement with 
the fact that the total mass of the K meson is considerably larger than that of the pion. 

Thus, formulae (4) and (5) with the a?s given in Table 1 suggest an interpretation 
of the main part of the rest mass: the mechanical mass of hot particles is an effect of 
their rotational motion in a three-dimensional abstract space whereby |a] is the “length 
of the angular momentum vector of this rotational motion. 


5. A screw model of hot particles 


The considerations of Sections 3 and 4 lead to similar conclusions: the attribute 
has something to do with a rotational motion in an abstract space. It is remarkable 
that these conclusions have been derived from quite heterogeneous arguments: in 


f 
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Section 3 we have taken account of the necessity of a strict conservation (of the third 
component) of the total angular isomomentum, whereas in Section 4 we have taken 
account of the increasing of |a| together with the mass of the corresponding particle, 
and particularly of the *magic^ ratios (3). 

However, there is also a serious discrepancy between the implications of Sections 3 
and 4. The considerations of Section 3 suggested the attribute to be a third component 
of an angular momentum, whereas those of Section 4 suggested that (the absolute 
value of) the attribute denotes rather the length of an angular momentum. The discre- 
pancy between these implications will be removed only if the third component is simul- 
taneously an extremum component of an angular momentum (maximum component 
for a particle and minimum component for an antiparticle). 

This suggests the following model of the particle: it is to be represented by a pair 
of points O and P in a three-dimensional euclidean space to be called isospace. The 
point O denotes the particle centre (analogously to the proton in the hydrogen atom) 
and the point P is able to rotate about the centre (analogously to the electron.) Both 
constituent points can possess spins of their own in isospace (analogously to the proton 
and electron in the atom). Both spins are simultaneously integer or half integer so 
that the total angular isomomentum of the hot particle is integer. The spin of the cen- 
tral point O will be interpreted as the isobaric (isotopic) spin of the hot particle while 
the angular momentum of the outer point Р (composed of its own spin and its orbital 
angular momentum with respect to О) will be connected with the attribute. 


In order to account for further details it will be necessary to introduce a privileged 
direction inside the particle. This direction may be joined organically to the particle 
centre by assuming the point Ото be not an ordinary point but a ‘directed point. (We 


denote the directed point О by O and the unit vector in its direction by о). A directed 
point in a three-dimensional space is a geometrical object with five degrees of freedom: 


three determining its position, and two its direction. The direction о is to play the role 
of “the electromagnetic third axis“ in isospace. The central point can possess a spin t 
oriented quite independently of the direction o but its component upon the privileged 
direction t.o is to play the role of a “third“ component tz. 


The outer point P is not a directed point but is subjected to two constraints: (i) it 
is bound to rotate at a given distance from the particle centre O, (ii) it rotates so that 


its momentum is parallel to the privileged direction o. 
If P is a spinless point, we get a simple classical model of the particle in isospace: 


the particle is composed of a directed central point O and of an outer point P rotating 


about O in a plan eperpendicular to the direction o in the sense of a right handed screw 
(Fig 1). For the antiparticle we assume a reversed situation: the point P is bound to 


rotate in a plane perpendicular to o in the sense of a left handed screw (fig. 2). In 


this way the projection a of the angular momentum of the point P upon the direction ó 
is identical with the length of the angular momentum vector |a| for particles and with 
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— |а| for antiparticles. This removes the discrepancy between the results of Sections 3 


and 4, 
Disregarding the spins of the points P i O our model possesses six degrees of free- 
dom (five describing the directed centre O and the sixth describing the position of the 


point P in a plane perpendicular to the electromagnetic axis 0). This model possesses 
a "size" given by the distance |OP| and a property of a right (or left) screw for particles 
(or antiparticles) respectively. By going over from a right-handed to a left-handed 


0 = 
(e Q 
= Г 
о РАНЕ 
С^ 
Fig. 1 Fig. 2 


system oÍ reference the roles of particles and antiparticles interchange. Thus, the 
connection between the operations of charge conjugation and inversion in isospace is 
quite analogous to that im the usual х-зрасе. 


Let us discuss a transition from a classical to a quantum model (disregarding still 
the spin of the point P). This transition, consisting in quantizing the orbital motion 
of P, is not unique. At first sight it seems that we have to quantize a plane rotator i.e. 


that we have to fix the position of Р in a plane perpendicular to the vector о. However, 
we may as well quantize the motion of P on a sphere and realize the supplementary 


condition (requiring the momentum of Р to be parallel to о) by choosing that linear 
combination of the 2|a| + 1 states (for which the square of the angular momentum 
is |а|(|а| + 1) for which the component of the angular momentum in the direction 


of o is an extremum (a maximum for a particle and a minimum for an antiparticle). 


Assuming a special frame of reference with the third axis oriented in the direction 0, the 
solution consistent with this supplementary condition is simply Yf, if the attribute of 
the particle is a. According to Formula (4), involving the square a(a -+ 1), we choose 
the latter alternative of realizing the supplementary condition in'a quantized version 
of the theory. 

If the point P possesses also a spin, the model is to be generalized so that (not 


its orbital but) its total angular momentum component is to assume a maximum 


(minimum) value in the direction o. It should be stressed that the supplementary con- 
dition concerns only the point P but not the point О. The spin of the point О (isospin) 
may be oriented quite arbitrarily. Therefore also the total momentum of the particle 


does not need to be oriented in the direction o. 
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6. A generalized postulate of a contact interaction 


In the foregoing section the privileged “electromagnetic” direction о was regarded 
as an intrinsic property of the particle. This is satisfactory since it deprives the isospace 
of anisotropic features. We may assume that the directions o for individual particles 
are distributed at random. This assumption will not violate the algebraic additivity 
of charge and attribute provided we introduce a generalized postulate of a contact 
interaction. This postulate may be stated as follows: 


Two or more particles do not interact directly one with another unless (i) their 
centres O and (ii) their directions OP coincide. 
The coincidence of the directed points O means that not only the points O but also 


their directions o coincide. In the traditional field theory the postulate of contact 
interaction could be secured by means of a set of Dirac delta functions and (upon 
integration) the coordinates of all individual particles or fields could be reduced to 
a single set of variables x, y, z (common for all fields or particles). Similarly now 
(in consequence of the generalized postulate of a contact interaction) the coordinates 


of all particles may be replaced by a common set of variables ифо ® denoting the 
range of variability of all (coincident) particle centres О, and two further variables Ф, 9 


denoting the range of variability of all coincident directions OP. We may assume further 
that neither the common position of the centres of the set of interacting particles nor 


their common direction o matter, i.e. cannot be distinguished from any other position 
or direction. Therefore, we are allowed to put (purely conventionally) the common 
centre Oat the origin of a Cartesian frame of reference and to fix (purely conventionally) 


the common vector o in the direction of the third axis. In this frame of reference the 
positions of the outer points P are described by means of their Cartesian coordinates 


Ujug ог by means of their spherical coordinates r, ф, 9. Since the direction OP is 
common for all interacting particles, we have a single set of variables p, 9. But the 
postulate of a contact interaction does not refer to the distance |ОР|, so that each particle 
may assume its individual r. On the other hand, since the physical objects situated 
at the points P are rotators, their individual r-values are not variable but constant. 
Therefore, the role of the variable r is different from the remaining variables. Disregar- 
ding this variable the field quantities are of the type y,, (p, 9) where а is an index 
of the isobaric spin (of the spin of the point O), and o is an index of the spin of the 
point P. In general Ya constitute а basis of a reducible representation of the group 
of rotations in isospace. Since the privileged axis has been identified (conventionally) 
with the third axis of a right-handed frame of reference, the field quantity must be 
of such a form that the angular momentum of the point P is parallel (antiparallel) 
to the third axis for each particle (antiparticle). 
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7. The meaning of the distance from the origin 


As pointed out in the foregoing section, the description of the particles in isospace 
may be reduced to an investigation of their dependence upon the spherical angles p, # 
and upon a distance r from the origin, ф, # are common variables for all interacting 
fields (particles), whereas the role of r is quite different. The values r for two or more 
interacting particles do not need to coincide, so that 7 is not a field coordinate in the 
usual meaning of this word. The question arises what is the physical meaning of the 
particle radius r in isospace? We shall fix its physical meaning by identifying it with 
the rest mass of the particle, In order measure r in conventional units of length, we 
have to introduce a constant / with dimension of length and assume 


Г ==. (6) 
Formula (6) shows that the variability of the point Р must be restricted to a surface 
of a sphere if our model is to describe a particle, i.e. a physical object with a definite 
rest mass. Thus, formula (6) is not only an explication of the physical meaning. of the 


distance from the origin, but may be regerded as a supplementary condition converting 
the free point P into a rotator. 


8. Field equations for integer isospins 


In order to guess plausible forms of field equations, let us start with an academic 
example of a particle being not a rotator but a free point P in isospace. Let us introduce 
the following abbreviations 


д ар О) 
=—i—, ---і--. (7) 
de д x, д1 


For a particle being a free point in both x-space and и-зрасе it is suggestive to assume 
a generalized Klein-Gordon equation 


(ри + E?) v (x, u) = 0 (8) 
This equation is separable into the usual Klein-Gordon equation and an eigenequation 
for the mass parameter 


(р? + т?) у (a) =0, (8-- m?) y (u) = 0. (8) 


These equations surely apply to the case of integer spins in both x-space and u-space. 
Since for an integer (half integer) spin of the point P also the isobaric spin (spin of 
the point O) is integer, (8) applies to the case of integer isobaric spin. If the particle 
is not a boson but a fermion, then there are two possibilities: we can start with (8) 
and, after separating it into (8), we can linearize the first equation (8’) while leaving 
the second as it stands: 


(iy, Р, +m) = 0, E (k? — m?) x = 0, (8") 
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or we can abandon (8) completely and assume a mixed equation 
Gy, p, + UP) y = 0. ©) 


where J is a constant with dimension of а length. (9) may be separated into 
(iy, Bie te mies (2 — m) = 0. (9') 


Thus, there is a single possibility of a mass eigenequation for bosons (the second 
equation (8’)) but two different possibilities of a mass eigenequation for fermions 
(either the second equation (8") or the second equation (9^). 

The mass operators appearing in (8) and (9) possess continuous spectra of eigen- 
values. This is connected with the academic assumption of a free point P in isospace. 
However, the point P is not a free point but a rotator. Taking account of the restriction 
(6), we have to modify the equations by supplementing them with suitable terms taking 


account of the reaction forces brought about by the constraints. The operator k? re- 
duces to 


hk? — E? = =; (10) 
Т 
where Ё is the radial momentum 
1 1 -- -— ili 
——[— — 11 
k sassa =| (11) 


and r= Іш. The operator 5? denotes the squared orbital angular momentum with 
eigenvalues а(а + 1) for a = 0,1... The mass eigenequations appearing іп (8”) and (8") 
or (9’) become | 


[8 — 12) — ту = 0, ог [KR — k) ту = 0, (12) 
whence we get, with the aid of (6) and (10), 


ты = a (a= Dior ты = T fa(at № а-0л,.. (08) 


The first formula (13) applies to bosons and fermions, the second only to fermions. 
In this way we have found an explanation of formula (4) for A and 2 as well as of 
formula (5) for л. 


9. Field equations for half integer isospins 


Since both the isospin and the spin of the outer point P are simultaneously integer 3 
or half integer, we have to linearize the mass eigenequation in order to describe particles 
with half integer isospins. A properly linearized second equation [8’) is 


(icok — m) y= ШІ = (14) | 
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where c are Pauli matrices for the point P and ø is the radial component of the vector с 


би 
o= —. (15) 


Е 
(14) applies to the academic example of a free point Р. For a rotator we have 


(oak - k) — m) y = 0. (16) 


Since o commutes with k but anticommutes with о k — ok, this equation may be 
written also in a more elegant form: 


і (с, 5 —m) =o. (17) 


This equation may be fused with the usual Dirac equation in to 


(va pat d es ZIP o (18) 


obviously applying to fermions. For bosons we may either iterate this equation which 


E 2 ШЕ 22 | v=o. (19) 


yields 


or use a mixed equation 


5 + — Tan ZIP = 0. (20) 
The eigenvalues of the operator 1 (с, c k] are Н “ағ a for a = JL : o 1 
2 г 2 Dek? 
whence, with the aid of (6), we get from (18) and (19) or (20) 
Т 1 "le - 1 T Уа 24459 Я 
„= (eri) or m= ydet) A CAT R GP US (21) 


where the first formula applies to fermions as well as to bosons, the second formula 
only to bosons. Formula (21) agrees with (4) for the baryons N, Е and with (5) for 
the meson К. However, the value / for bosons is about twice as large as that for baryons. 


10. The hypothesis of isotrino 


In the foregoing sections we were concerned only with the hot particles listed in 
Table 1 and disregarded the leptons. Under this restriction a conservation of a total 
angular isomomentum has been established not only for strong but also for weak 
interactions. Strong interactions are characterized by a conservation of both isospin 

£ 
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and attribute (separately). This means that strong interactions are characterized by an 
interaction term in form of a product of two invariants: one of them involving merely 
theisospins of the fields in question (i.e., the spins of the central points О of all particles), 
the other involving merely the angular momente of the outer points P. On the other 
hand, weak interactions between hot particles do not conserve the isospin and the 
attribute separately but do conserve the total angular isomomentum. This means that 
there is an exchange of isomomentum between the central point Oand the outer points 
of the particles in question. Such exchange can be described by means of interaction 
terms involving invariant products of two vectors: the isospin vector and the angular 
momentum vector of the outer points P. In this way the invariance under rotations 
in isospace is secured even for weak interactions between hot particles. 

However, the above described device fails for interactions involving leptons, e.g. 


for the pion decay z — и + v whereby a total angular momentum |t| = 1 disappears 
(the particles to the right possessing neither an isospin nor an attribute). This difficulty 
resembles the situation after the discovery of the beta decay. There was also an apparent 
violation of conservation laws but this difficulty has been overcome by the neutrino 
hypothesis of Pauli. Similarly let us assume the existence of a new particle, the isotrino, 
possessing neither energy-momentum nor angular momentum but possessing an 
isomomentum. Thus, the isotrino exists only in isospace, i.e., is described by a field 
quantity independent of the x-coordinates. If it is so, we have three types of particles (i) 
hot particles existing simultaneously in both spaces (x-space and isospace), (ii) leptons 
existing only in x-space, and (iii) isotrinos existing only in isospace. In this way the 
roles of both spaces become more symmetric. 


We may assume that isotrino possesses an кошо a vanishing attribute. 
The decay of the pion should involve the emission of a pair of isotrinos. The weak 


ң і : А г 
interactions between hot particles characterized by Да = > may be also regarded аз 


involving isotrino: the exchange of isomomentum с the central point О and 


the outer point P may be regarded as an exchange of a virtualisotrino between Qand P 
instead of a direct interaction. This would account for similar orders of magnitude of 
both types of weak interactions: with and without a well defined Да and At. 
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Г. Райский, О винтовой модели частицы в изопространстве 

Рассматривается в изо-пространстве модель эгорячих” частиц, которая позва- 
ляет на физическую интерпретацию заряда, изобарного спина, странности и мас- 
сы покоя как следствия оборотного движения в изо-пространстве. Эта модель 
имеет некоторое сходство с атомом водорода: частица имеет ядро со спином, 
обращающийся вокруг него второй пункт имеющий спин, орбитальный момент 
количества движения и определенную величину. Частица (античастица) имеет . 


- 


au sd L =. 
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тоже правовинтовой (левовинтовой) винт. Спин ядра интерпретируется как спин 
изобарический. „Странность” заменяется более физическим атрибутом будущим 
зкстремальной составляющей момента количества движения наружной точ- 
ки. Этот момент количества движения, вместе с радиусом частицы, определяет 
| массу покоя. Привилегируемая „третья ось” не есть признаком изо-пространства, 
но самой частицы. Алгебраическая аддитивность странности обезпечона благо- 
| даря предположению контактного действия. Полный изо-момент количества 
движения есть величиной постоянной как для сильных так и для слабых воз- 
| действий между „горячими” частицами. Существование новой „частицы — изо- | - 
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LABORATORY EQUIPMENT AND TECHNIQUES 


AN AUTOMATIC TIMING AND RECORDING CIRCUIT FOR USE WITH 
A NEUTRON CRYSTAL SPECTROMETER 


By Denis О’Соммов AND Konrap BLINOWSKI 


Institute of Nuclear Research, Polish Academy of Sciences, Warsaw. 
(Received February 28, 1957) 


Apparatus is described for the automatic control and recording of the measurements 
оп a crystal neutron spectrometer. The angular position, time of measurement and intensity 
of the prumary and scattered neutron beams are registered and recorded on ап electric 
desk adding machine. The circuits ensure the most economic use of the measurement 
time available by relating the time duration and angular frequency of measurement points 
to the scattered beam intensity. 


In making measurements with a neutron spectrometer or a neutron diffracto- 
meter it is necessary to record the neutron intensity of the diffracted beam as a function 
of the angle of diffraction. This may be accomplished in two ways: (a) the neutron 
counter can be rotated about crystal at constant speed and the counting rate of the 
counter recorded by a counting rate meter driving a pen recorder (b) measurement 
is made on a scaler with the counter stationary, the count on the scaler is recorded 
and the counter is then moved on to the next measurement position. The circuit 
described below is designed for recording by the second of the above methods. 

In using this second method it is necessary to preset both the time measurement 
interval and the angular interval between readings. If two quantities are fixed time will 
be wasted counting high intensity peaks for an unnecessarily long time, or taking too 
frequent readings in angular regions of the diffraction pattern which are of low intensity 
and of small interest (for example between the peaks of a diffraction pattern). In the 
present design the measurement time can vary within wide limits which are preset. 
High intensity peaks are automatically counted for a shorter time than lower intensity 
points, the time being varied to keep the statistical accuracy of the readings more or 
less constant. Similarly the angular interval between readings is automatically reduced 
at regions of high counting intensity in the diffraction pattern. š 
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The basic feature of the circuit is the use of a standard electric desk adding machi- 
ne to print the values of all the measured quantities in decimal form. The following 
quantities are measured and recorded: 


l. recorded counts by neutron counter 
2. time taken for measurement 

3. count rate of monitor neutron counter 
4. angular setting of spectrometer arm. 


Impulses in the two neutron counters are registered by hard valve three decade 
scalers. The output 1000-th impulses are counted by 1000-th registers (see Fig. 1). 
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Fig. 1. Block schematic diagram of automatic recorder. 4 — Angle Divider, B — Angle interval limit, 

С — Angle Register, Е — Convertor, F — Counter, С — Gate, H — Printer, J — 1000 Register, 

J — Sample Drive, K — Scaler A, L — Scaler B, M — Sequence control, N — Spect. drive (arm), 
O — Time interval limit, P — Time Register. 


Fig. 2. Impulse stretching circuit. PA — telephone relay type B, resistance 2000 ohms. WPA — сой 
of uniselector relay. Condenser connected to grid circuit by contact РА charges through resistance 
maintaining positive potential on valve grid after positive impulse from scaler has already decayed. 


Telephone stepping relays (uniselectors) are used as registers. The position of the wiper 
arms of the stepping relays depends on the number of impulses received by the relays. 
The contact banks of the stepping relays are connected in appropriate order to ten 
electromagnets mounted on the keyboard of the adding machine. When the wiper 
arms are earthed one of the ten magnets is energized, according to the position of the 
wiper arms and hence corresponding to the number of impulses stored in the relay. 
Each 1000 register has a capacity of 99 impulses. A valve circuit *stretches* the output 
1000 th impulse from the scalers to a time duration sufficient to operate the uniselector 
(see Fig. 2) A typical example of a register circuit, in this case the time register, is 
given in Fig. 3. 

Timing impulses are supplied to the circuit from a pendulum clock and are stored 
in a time register also based on a stepping relay. These impulses, one per minute, also 
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control the starting and stopping of the scalers via the “sequence control“ circuit which 
is described below. 

The angular position of the spectrometer arm is read at the spectrometer itself 
by the aid of a divided and figured glass circle with a scale division for every minute 
of arc. For automatic recording a photoelectric device detects the passage of the scale 
divisions across the fiducial mark and provides an electrical impulse for every minute 
of arc rotation. These impulses are fed to an angle register consisting of 4 stepping 
relays and having a total capacity of 9.999 minutes. These relays have their appropriate 
contacts connected to the adding machine electromagnets and in addition switch on 


o«60V 
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Fig. 3. Time register circuit. ИМС — coil of uniselector. ВІ, B2, and B3 — contact banks of uniselector, 

dm — interrupter springs of WMC. 0,1,2..., 9, + — electromagnets operating keys of adding machine. 

С — minute impulse relay contact. ЯМХ — reset relay. When the time measurement is printed contacts 

1, 2, and 3 are earthed in that order by the sequence control circuit. Energizing the + electromagnet 

causes the number registered on the uniselector to be printed. Contact 4 is normally connected to + 60 V, 
when momentarily disconnected uniselector WMC is reset to zero. 


telephone indicator lamps to give a visual indication of the angle reading. 40 lamps are 
connected in 4 decades of lamps labelled 0 to 9 giving a decimal reading of angle. 
The stepping relays can be set to correspond to the position of the spectrometer arm 
before starting a series of measurements and will thereafter keep in step with the 
movement of the arm. т & 

The operation of the whole circuit is controlled by the circuit labelled ,,sequence 
control“ in the block diagram. This circuit is based on a stepping relay. On receipt 
of the first minute impulse the sequence control circuit operis the gates A and В, 
the gate between the clock, and the time register. Scaler A registers the impulses from 
the monitor counter. On the next minute gate Æ is closed so that scaler A and its 1000-th 
register give a reading in impulses per minute of the intensity at the monitor counter. 
The latter is so placed as to give a sufficiently large counting rate to ensure good 
statistical accuracy for a 1-minute count. Scaler B continues to register impulses until 
stopped by a time interval limit circuit (see Fig. 4). This circuit is based on an idea 
suggested by Taylor (1950). Timing impulses and 1000-th impulses from scaler B are 
fed together to an adding circuit based on a stepping relay. When the sum in the 
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adding circuit attains a certain value an auxiliary circuit is prepared which switches 
off the scaler on receipt of the next minute impulse (or immediately if the last impulse 
to the adding circuit was a minute impulse). The effect of this circuit is to automatically 
adjust the measurement time to the counting rate. If 1000 n impulses are counted 


in m minutes then n and m fulfil the condition 


n+m=N 


Fig. 4. Time interval control circuit. WO — coil of uniselector. B1 and B2 — contact banks of uniselector. 

C — contact of minute impulse relay. PB — contact of 1000th impulse relay. WOX — reset relay, mode 

of action as WMX in Fig. 3. KP — prepare finish relay. SB — start scaler B relay. P — switch selecting N. 

On earthing contact 1 momentarily, relay SB operates and hold: by itd own contact. C and PB impulses 

cause WO to rotate and earth one coil of KP which is held by its own contact and the second coil. The 
next C impulse then causes SB and KP to return to their inactive condition. 


when N can be chosen to suit the operating conditions. With № = 10 the time m for 
various counting rates n/m is as follows: 


n/m m statistical error % 
9 1 1.05 
5 2 1.10 
Я. 5 1.40 
0.5 10 3.30 


It is seen that while the counting rate and time vary within wide limits the statistical 
error is held relatively constant. 

When scaler В is stopped by the time interval limit circuit the sequence control 
unit starts the printing cycle (the angle reading is printed first). The sequence control 
unit is based on a stepping relay which, in rotating, earths in turn the wiper arms of 
the angle register stepping relays. Thus when the wiper arm on the 10,000 ths relay 
is earthed one of the electromagnets on the printing machine is energized corresponding 
to the number stored on that relay. Then in turn the number of 1000s, 100s 10s, and 
units is registered by the machine and then printed. When the angle has been registered 
the sequence control unit sends out a signal to the spectrometer arm drive unit, the 
operation of which is described later. 

The printer then proceeds to print the time reading from the time register and 
the readings on the two scalers and their 1000-th registers. After printing the time 
reading the reading on the scaler A 1000-th register is registered by the printer. Next ке 
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the 100s, 10s, and units readings of the decades must be registered by the machine 
before the whole number is printed. This is carried out by the aid of a special conversion 
unit (see Fig. 5). Although the scaler used is a decade scaler each decade is in fact 
based on four scale-of-two circuits using the 1-2-2-4 system. Each decade has four 
output leads connected to the four scale-of-two circuits and these four outputs are 


Fig, 5. ТІ... ТА contacts of relays T1... T4. Relays, T1, T2, T3, and ТА are operated by the scales of 

two in a decade scaler based on the 1-2-2-4 system to which they are connected’ by circuits shown in 

5 (b). Contacts DA1, DA2, DA3, DB1, DB2, DB3 connect in appropriate order the valve circuit to the 

scale of two valves in six separate decades. Four circuits 5 (b) are used for the relays Т1... T4, the latter 

thus reflect the state of the 1-2-2—4 scale of two circuits in the decades. Earthing contact 1 operates one 
of the ten electromagnets according to the state of the decades. 


connected to a valve and relay circuit which converts the 1-2-2-4 system to a decimal 
output. Only one converter circuit is used which is connected to the decades in appro- 
priate order by means of gating circuits. 

When the printing cycle has been completed the scalers, 1000-th registers, time 
register, and time interval limit circuit are reset to zero. Before the circuit is restored 
completely to its initial condition it waits for the spectrometer arm to complete its 
movement. The angular movement of the arm between readings is controlled by ап 
“angle interval limit“ circuit based on a similar principle to that time interval limit 
circuit. ; ⁄ s 

If 40 is the angular shift then it must fulfil the following condition 


na 40 : 
PEL Dep | 


where n/m is the intensity of the neutron beam at the last measurement and z and Ө 
are constants which can be preset within a wide range of values. In general if the neutron 
intensity of the diffracted beam increases then 40 will decrease, i.e., measurements 
will be made at smaller angular intervals. This feature should be useful when making 
powder diffraction pattern measurements. By making @ large and x small 40 can 
be made more or less independent of n/m if it is so desired. When the arm begins to 
move, after receipt of the arm drive start signal (see above), minute of arc impulses 
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from the angle measurement system are fed through an angle divider circuit to the 
angle interval limit circuit. The angle divider circuit, based on a stepping relay, gives 
output impulse for every y input pulses. y can be set equal to 1, 3 or 9. The output 
impulses /0/y are mixed in the angle interval limit circuit with n/m impulses stored 
there during the counting cycle. The latter are derived by feeding in 1000-th impulses 
from scaler B during the interval of operation of scaler A, i-e., 1 minute. When the 
sum (40/%) + (n/m) reaches the present value Ө the angle interval limit circuit stops 
the spectrometer drive. This stop signal is simultaneously fed to the sequence control 
circuit which is then restored to its initial condition. On receipt of the next clock impulse 
the whole cycle is repeated. 

For use in transmission measurements a sample drive mechanism is mounted 
on the spectrometer arm. This mechanism is controlled by the automatic controller. 
When in use, each ,,start arm drive“ signal is fed to the sample drive motor and simul- 
taneously throught a relay scale of two circuit to the spectrometer arm drive circuit. 
Thus the sample is moved twice for every movement of the arm, i.e., into and out of 
the beam. 

An additional circuit, not shown in the diagram, enables the automatic reduction 
of 40, by reduction of y, at up to 4 angular positions which can be preset through 
the whole angular range. 

The authors wish to thank Professor Bronislaw Buras for his interest and en- 
couragement afforded during the course of this work. Thanks are also due to Stefan Bo- 
chenek, Ryszard Kula, and Aleksandra Skarzyüska for assistance in building the 
apparatus. 


KPATKOE СОДЕРЖАНИЕ 


Jl. О’Коннор, K. Блиновский: Автоматическое устройство для управления 
неитронным кристаллическим спектрометром. 


Описанное устройство предназначено для регистрации измерений и автома- 
тического управления нейтронным спектрометром. Регистрируется угловое поло- 
жение, время измерения, интенсивность первоначальных и рассеянных нейтро- 
нов. Для этого использовано типовой электрический сумматор. Устройство дает 
возможность использовать наиболее экономически время измерений по той при- 
чине, что время и частота измерений зависят от интенсивности рассеянных ней- 
тронов. 
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It is well known that the standard molecular orbital method is not very successful 
when used to calculate the excitation energies of molecules. Moffitt (1951) and Pariser 
(1953) have pointed out that a major source of error in these calculations is the in- 
accuracy of the energies of asymptotic dissociation products of molecular states, and 
the above authors have also given the ways of eliminating this error. Pariser considered 
two isolated carbon atoms in their sp? valence states and the C+ and C- ions. In the 
first case the energy of two 2pz-electrons is given by 2W,; in the second case, when 
both electrons are attached to the same atom, the energy of the ion pair is 2/7, + 
+ (аа | aa), where (аа | aa) is the one-centre interaction integral 


(aa | aa) = Í f ф 0) 92 Qn) pa) 9, (2) dr, des (1) 


in which ф, denotes the 2pz atomic orbital of the a-th atom. The difference of energy 
between the ion pair and the neutral pair is then simply (аа | аа). Experimentally 
however this difference is given by the ionization potential J of a neutral carbon atom 
in its valence state minus the electron affinity A of the same atom. ‘Thus we should have 


(аа | aa) = I, — A,. (2) 


Using experimental values of the ionization potential and of the electron affinity 
of a carbon atom in its sp? valence state we get (aa | aa) = 11.54 — 0.46 = 11.08 eV, 
whereas the theoretical value of this integral is 16.93 eV. | 

An essential point of the successful theory developed by Pariser and Parr (1953a, b) 
is the use of the semi-empirical values given by (2) for the (аа | aa) integrals instead | 
of the theoretical ones which can be obtained from formula (1). However, it has been 


(299) 
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300 š W. Kolos 


shown (Kolos 1957a, b) that the results of the molecular orbital method can be also 
improved by considering the effect of the correlation of electrons; instead of the usual 


two-electron wave functions @,, = %,(1) v,(2) the functions | 
T єў ape Vn (1) Yn (2) Т1 (3) 


have been proposed. 
If wave functions of the form (3) are used to calculate the (аа | аа) integral, 


then this integral is given by 
f #00202) тара) ра) dni dz 
Г] воде) фа) ға (2) ат, dz, 


The values of the (aa | аа) integrals calculated from (1), (2) and (4) are set out 


in Table I. 


(4) 


(аа | аа) = 


Table I 


Values of the One-Center Interaction Integrals 


10.533; 11.084 


12.278 
14.525 


1 Z denotes the effective charge in the Slater-type atomic orbital. 

2 The figures in brackets refer to the formulas used in the calculation. 

3 Pariser and Parr 1953b. 

4 Parr and Pariser 1955; in this paper new values of the ionization potential and electron affinity 
were used. 


5 Fumi and Parr 1953. 


We see that the wave functions (3) give results in agreement with those which can 
be obtained from experimental data and which were used in the semi-empirical. theory 
of Pariser and Parr. We must, however, be very careful when drawing conclusions from 
this agreement. Formula (2) was obtained from the difference of energy between the 
ion pair and the pair of neutral atoms, and it does not hold it two-electrons functions (3) 
are used for the ion. It has been also shown by Pauncz (1954) that the procedure of 
separating the problems of atomic from those of molecular structure is not free from 
objections. Pauncz has calculated the ground state energy of фе H, molecule; in 
the ionic terms of the molecular orbital wave function he has used the function giving 
the correct energy of the H- ion. However, in this way he obtained a higher energy 
for the H, molecule than that resulting from the use of simpler functions for the ionic 
terms. Nevertheless, it seems that the results presented above throw some new light 
on Pariser and Parr's theory, and show the physical meaning of using the semi- 
empirical values of the one centre interaction integrals in the molecular orbital method. 


BEF zm 3 
The One-Centre Interaction Integrals _ 
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